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ABSTRACT 

Clusters of galaxies have long been used as laboratories for the study of galaxy evolution, but despite 
intense, recent interest in feedback between AGNs and their hosts, the impact of environment on these 
relationships remains poorly constrained. We present results from a study of AGNs and their host 
galaxies found in low-redshift galaxy clusters. We fit model spectral energy distributions (SEDs) to the 
combined visible and mid-infrared (MIR) photometry of cluster members and use these model SEDs 
to determine stellar masses and star- formation rates (SFRs). We identify two populations of AGNs, 
the first based on their X-ray luminosities (X-ray AGNs) and the second based on the presence of a 
significant AGN component in their model SEDs (IR AGNs). We find that the two AGN populations 
are nearly disjoint; only 8 out of 44 AGNs are identified with both techniques. We further find that 
IR AGNs are hosted by galaxies with similar masses and SFRs but higher specific SFRs (sSFRs) than 
X-ray AGN hosts. The relationship between AGN accretion and host star-formation in cluster AGN 
hosts shows no significant difference compared to the relationship between field AGNs and their hosts. 
The projected radial distributions of both AGN populations are consistent with the distribution of 
other cluster members. We argue that the apparent dichotomy between X-ray and IR AGNs can be 
understood as a combination of differing extinction due to cold gas in the host galaxies of the two 
classes of AGNs and the presence of weak star-formation in X-ray AGN hosts. 
Subject headings: galaxies: active, galaxies:clusters:general, infrared radiation, X-rays 



> 

(N 
00 

o 



X 



1. INTRODUCTION 

Galaxy formation and evolution has long been a sub- 
ject of considerable interest, with early work dedi- 
cated to explori ng the physical processes responsible for 
star- formation (I Whipple! 1 1 9461). expla ining the genesis 
of the Milky Way ( Eggen et al.l Il962t ) , and examining 
the e volution of galaxies in clusters (jSpitzer fe B aadc 
1951). Models for the evolution of galaxies in clus- 
ters gained strong observational constraints with the 
discovery of a n apparent ev olutionary sequence among 
local clusters (jOemlerl 119741 ) . The discovery that the 
fraction of blue, spiral galaxies in relaxed galaxy clus- 
ters in creases from z = to z « 0.4 quickly fol- 
lowed ([Butcher fc Oemlerl fl978l U98l . The dearth of 
spiral galaxies in the high-density regions at the cen- 
ters of galaxy clusters is known as the morphology- 
density relation HDresslerlll980t IPostman fc GeUerl ll984l; 
iDressler et all Il997t iPostman et al.l 120051 ). This rela- 
tion places additional, strong constraints on evolutionary 
models for cluster galaxies. That star-forming galaxies 
are also rare in the centers of clus ters had be e n pre - 
viously suggested by the results of lOsterbrockl (1960) 
and was subsequently ob served in other work (Gislcr 
119781 : IDressler et al.l 119851 ) . The impact of environment 
on the frequency and intensity of star-formation at a 
wide variety of dens ity scales has been measured us- 
ing numerous visible ([Abraham et al. l 1996: Balog bTet al.l 
[R^lKauffmann et al.ll2004t iPoggianti et al.ll2006ll2008fc 



von der Linden et al. I2010D and m id-infrared (MIR; 
Saintonge et al.ll2008t iBai et al.ll2009f ) diagnostics. Star- 
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forming galaxies are consistently found to be more com- 
mon and to have higher star-formation rates (SFRs) 
in lower density environments and at high e r re dshift 
(jKauffmann et al.|[200l IPoggianti et al.ll2006l I2008I ). 

The observed trends in star-formation with environ- 
ment are usually attributed to variations in the sizes of 
gas reservoirs, either the existing cold gas or the hot gas 
that can cool to replenish the cold gas as it is consumed. 
Given that AGNs also consume cold gas to fuel their 
luminosity, similar patterns might be expected among 
AGNs. Indeed, recent work reveals strong dependen- 
cies of the l uminosities and types of AGNs on envi- 
ronment (e.g. jK auffma nn et al.ll2004l: iPopesso fc B viand 
120061: iConstantin et alJl2008t iMontero-Dorta et al.ll2009t ) 
for AGNs selected via visible-wavelength emission-line 
diagnostics. Von der Linden et al. (2010) find fewer 
"weak AGNs" (primarily LINERS) among red sequence 
galaxies near the centers of clusters compared to the 
field, but they find no correspondin g dependence amon, 
blue g alaxies. Intriguingly, while IMontero-Dorta 
(2009) independently report a decline in the fraction of 
low-luminosity AGNs toward the centers of low-redshift 
clusters, they find an increase in the fraction of LIN- 
ERs in higher density envi ronments. The differ e nce is 
likely a result of evolution. IMontero-Dorta et al.l (|2009l ) 
found qualitatively different behavior between their main 
z ~ 1 sample and the result produced when they applied 
their analysis to SDSS clusters. These results indicate 
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that the variation of galaxy properties with local environ- 
ment may influence the types of AGNs observed and that 
evolution in the relationship between some AGN classes 
and their host galaxies is important. Understanding the 
environmental mechanism that transforms star-forming 
galaxies into passive galaxies in clusters may help relate 
gas reservoirs in cluster galaxies to galaxy evolution as 
well as to AGN feeding and feedback. 

Several mechanisms to cause the transformation 
from star-forming to passive galaxies have been 
proposed. These includ e ram - pressure str i pping 
of cold gas dGunn fc Gottl fl972t lOuilis et all f200a 
Roediger fc Henslerl l2005l)~ ~ strangulation (lLarson et all 
19801: iBalogh et al l 12000 : Kawata fc Mulchaevi 120081 

McCarthy et al. 20081) and galax y harassment 

(jMoore et all 119961 11991 ILake et all 119981 ). Each 
mechanism operates on a different characteristic 
timescale and has its greatest impact on galaxies of 
different masses and at different radii. In principle, 
the transition of galaxy populations from star-forming 
to passive as a function of environment can probe 
the relative importance of these processes. However, 
such appr oaches su ff er fro m practical difficulties. For 
example, iBai et all {2009) argue that the similarity 
of the 24/xm luminosity functions observed in galaxy 
clusters and in the field suggests that the transition from 
star-formation to quiescence must be rapid, which im- 
plies that ram pressure stripping must be the dominant 
mechanism. Von der Linden et al. (2010), by contrast, 
find a significant trend of increasing star-formation with 
radius up to 5i?2oo from cluster centers. They conclude 
that preprocessing at the group scale is important, 
which is inconsistent with ram pressure s tripping as 
the d river of the SFR-density relation. iPatel et all 
(2009) find a similar trend for increasing average SFR 
with decreasing local density down to group-scale 
densities (T, ga i w 1.0 Mpc~ 2 ) near RX J0152. 7-1357 
(z = 0.83). The importance of preprocessing in group- 
scale environments reported by these authors suggests 
that strangulation rather than ram pressure stripping 
drives the SFR-density relation. The st arkly different 
conclusions reached by IBai et all (120091) comp ared to 
IPatel et all (|2009l) and Ivon der Linden et all (|2010l ). 
despite their common use of star-forming galaxies to 
examine the influence of environment, highlight the 
difficulties inherent in such studies. 

Attempts to distinguish between various environmen- 
tal processes become still more difficult with cluster sam- 
ples that span a wide range in redshifts. The epoch 
of cl uster assembly (0 < z < 1.5, e.g. iBerrier et all 
2009) coincides with the epoch of rapidly declining star 
forma tion (e.g. iMadau et all 119981 : iHopkins fc Beacoml 
20061) and AGN activ i ty Ce.g. iShayer et all 119961 
Bovle fc Terk^vichl [l99l iShankar et all 120091 ) . which 



and star-tormm g gala xies tFoggian ti et al.l I^UUO , IZUU8t 
iSaintonge et~aLl 120081: lHaines et all 120091 ) identified us- 
ing a variety of methods. These newer results have also 
examined cluster members confirmed from spectroscopic 
redshifts rather than relying solely on statistical excesses 
in cluster fields, which permits more detailed study of the 
relationships between galaxies and their parent clusters. 

The wide variety of AGN selection techniques em- 
ployed in more recent studies represents an impor- 
tant step forward in understanding the dependence of 
AGNs on environment. Several recent papers have 
used X-rays to study the frequency and distr ibution 
of AGNs in galaxy clusters (IMartini et al.l 120061 hence- 
forth M06; Martini et all l2007r ' 
lArnold et a l H2009l: lHart et al.lT2 009) 
with reds hift (Ea stman et alll2007t 



star-forming galaxies (iPoggianti et al.l 120061 120081 : 



makes it difficult to disentangle rapid environmental ef- 
fects from the global reduction i n the amount of avail- 
able cold gas. iDressler fc Gunnl (|1983l) found early ev- 
idence for an increase in AGN activity with redshift, 
and the Butcher-Oemler effect had already provided ev- 
idence for a corresponding increase in SFRs. In the 
last decade, the proliferation of observations of high- 
redshift galaxy clusters at X-ray, visible and infrared 
wavelengths h as yielded similar trends in the fraction 
of both AGNs (Ea stman et al.ll2007t IMartini et alll2009D 
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and their evolutio n 
Martini et aT1l2009l ). 



IMartini et "all (2009) found that the AGN fraction among 
cluster members increases with decreasing local density 
and increases dramatically (Jagn oc (1 + z ) 5 - 3±1 - 7 ) with 
redshift. They also found that X-ray identification pro- 
duces a much larger AGN sample than visible- wavelength 
emission line diagnostics: only 4 of the 35 X-ray sources 
identified as AGNs by M06 would be classified as AGNs 
from their visible-wavelength emission lines. Similar re- 
sults have been found when comparing radio, X-ray and 
mid-IR AGN selec tion techniques for field AGNs (e.g. 
IHickox et al.ll2009l) . 

The different AGN selection techniques identify dif- 
ferent AGN populations and suffer from distinctive se- 
lection biases. Both X-ray and visible-wavelength tech- 
niques can miss AGNs due to absorption, either in the 
host galaxy or in the AGN itself; however, X-ray selection 
can find lower luminosity AGNs and AGNs behind larger 
absorbing columns compared to emission line selection. 
Mid-infrared selection techniques suffer from relatively 
poor angular resolution, so they are mainly sensitive to 
AGNs that outshine their host galaxies in the band(s) 
used to perform the AGN selection. The X-ray and visi- 
ble techniques can also be contaminated by emission from 
the host galaxy. While the identification of X-ray sources 
with Lx > 10 42 erg s _1 as AGNs is unambiguous, X-ray 
luminosities in the 10 40 -10 42 erg s _1 range can be pro- 
duced by low-mass X-ray binaries (LMXBs), high-mass 
X-ray binaries (HMXBs) , and thermal emission from hot 
gas. Both visible-wavelength and MIR indicators are 
subject to contamination from young stars, which pro- 
duce emission lines and heat dust near star-forming re- 
gions until it emits in the MIR. Even the interpretation of 
the well-established Baldwin-Phillips- Terlevich diagram 
(jBaldwin et al.|[l98l can be controversial in the transi- 
tion region between star-forming galaxies and AGNs. 

These difficulties motivate the use of multiple tech- 
niques to obtain a complete census of AGN and to cor- 
rectly identify potential imposters. In this paper, we 
extend the work of Martini et al. (2006, 2007) by supple- 
menting their X-ray imaging and visible- wavelength pho- 
tometry with MIR observations from the Spitzer Space 
Telescope. We use these data to select AGNs indepen- 
dent of their X-ray emission. We also measure the prop- 
erties of AGN host galaxies by fitting their visible to 
MIR spectral energy distributions (SEDs). We discuss 
our visible and MIR data reduction and photometry in 
Section[2j Section[3]details our techniques for identifying 
AGNs and measuring galaxy properties, and we describe 
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the results in Section @] We discuss the implications for 
the relationship between AGNs and their host galaxies 
in Section Throughout this paper we use the WMAP 
5-year cosmology — a AC DM universe with f 2 m = 0.26, 
Q A = 0.74 and h = 0.72 (|Dunklev et al.ll2009| ). 

2. OBSERVATIONS & DATA REDUCTION 

We obtained MIR observations with the Spitzer Space 
Telescope of the X-ray sources identified as members of 
8 low-redshift galaxy clusters by M06. The initial reduc- 
tion of the Spitzer imaging is described in Section 12.11 
Visible wavelength photometry of these clusters were ob- 
tained at the 2.5m du Pont telescope at Las Campanas 
by M06. We provide a brief summary of these data in 
Section 12.21 further details are provided by M06. We 
then discuss the corrections for Galactic extinction and 
for instrumental effects in Section l2~3l 



2.1. Spitzer Reduction 

We obtained mid-infrared (MIR) observations from 
the Spitzer Space Telescope using t he IRAC (X e ff = 
3.6 , 4.5,5.8,8.0 u?7i:lFazio et al.| [200l and MIPS (A e// = 
24; iRieke et al.l 20041) instruments from Spitzer program 
50096 (P.I. Martini). Observations were carried out be- 
tween 2008 November 1 and 2009 April 22. Spitzer point- 
ings were chosen to image the X-ray point sources in 8 
low-redshift galaxy clusters identified by M06. We sup- 
plemented these observations with data from the Spitzer 
archive for Abell 1689 and AC 114. 

Spitzer's cryogen ran out before the MIPS observations 
of three clusters (Abell 644, Abell 1689 and MS 1008.1- 
1224) were carried out. In one of these clusters (Abell 
1689) we extended our coverage to 24/im using observa- 
tions from the Spitzer archive, leaving two clusters with 
no usable MIPS observations. The Astronomical Obser- 
vation Request (AOR) numbers used to construct the 
MIR mosaic images of each cluster are listed in Table 
1, along with the corresponding 3er observed- frame lu- 
minosity limits at both 8 and 24/im. These limits are 
approximate because the image depth varies across the 
mosaics due the changing number of overlapping point- 
ings. Quoted limits correspond to areas with "full cover- 
age" but without overlap from adjacent pointings. 

The raw Spitzer data are reduced by an automated 
pipeline before they are delivered to the user, but arti- 
facts inevitably remain in the calibrated (BCD) images. 
Preliminary artifact mitigation for the IRAC images was 
performed using the IRAC artifact mitigation tool by 
Sean CarejQ- We inspected each corrected image after 
this step and determined whether the image was imme- 
diately usable, if additional corrections were required, 
or if it simply had too many remaining artifacts to be 
reliably corrected. The latter class primarily included 
images with extremely bright stars that caused artifacts 
too severe to be corrected. Where appropriate, addi- 
tional corrections were applied using the muxstrip^3 and 
jailbaiQ correctors by Jason Surace and the column pull- 

1 http://spider.ipac.caltech.edu/staff/carey/irac_artifacts/ 

2 http: / / ssc.spitzer.caltech.edu/dataanalysistools / tools / contributed/ 
irac/ automuxstripe / 

3 http: / / ssc.spitzer.caltech.edu / dataanalysistools / tools / contributed / 
irac/jailbar/ 



down correctoiQ by Leonidas Moustakas. Artifacts in the 
MIPS images were removed by applying a Hatfield cor- 
rection algorithm packaged with the Spitzer mosaic soft- 
ware, (MOPEX0), as described on the Spitzer Science 
Center (SSC) websitfl 

Mosaic images for both IRAC and MIPS were 
constructed from the artifact-corrected images using 
MOPEX. Aperture photometry was extracted from the 
resulting mosaics using the apphot package in IRAF. We 
converted the measured fluxes to magnitudes in the Vega 
system after the photometric corrections described in 
Section 12.31 had been applied. All magnitudes quoted 
in this work, both visible and MIR, are calculated with 
respect to the Vega standard. The photometric apertures 
used by apphot were chosen to enclose a region of approx- 
imately 10 kpc projected radius at the redshift of each 
cluster. These large apertures yielded reduced S/N, but 
most cluster members were sufficiently bright that the 
uncertainties on the measured fluxes were dominated by 
systematic errors (5%) in the zero-point calibration, ex- 
cept at 24/im. The use of large photometric apertures 
also allowed galaxies to be treated as point sources for 
the purpose of computing aperture corrections, as rec- 
ommended by the SSC. A smaller aperture could im- 
prove the S /N, but this gain would be outweighed by the 
systematic uncertainty introduced by the aperture cor- 
rections for the resulting flux measurements, as aperture 
corrections for IRAC extended sources remain highly un- 
certain (IRAC Instrument Handboolfl) . 

2.2. Visible Photometry 

All 8 clusters in our sample have B-, V- and i?-band 
imaging, and 4 of the 8 have /-band imaging. We ex- 
tracted separate source catalogs for e ach of these bands 
using Source Extractor (SExtractor, iBertin fc Arnouti 
1996) and merged the catalogs using the i?-band image 
as the reference image for astrometry and total (Kron) 
magnitudes. We correct from aperture to total magni- 
tudes without altering the colors from the aperture pho- 
tometry by applying the i?-band aperture corrections to 
all bands, 



m Kron = m, Ap - (R/ 



R 



Ki 







(1) 



where rriAp and mKron are the aperture and Kron-likc 
magnitudes, respectively, for the band being corrected. 
Rather than taking the published photometry from M06, 
we used the redshift-dependent apertures assigned to 
each cluster as described in Section 12.11 This maintains 
consistency with our IRAC photometry and results in rel- 
atively small aperture corrections, typically ~ 0.1 mag. 

SExtractor returns i?-band positions that are good to 
within a fraction of an arcsecond. However, the posi- 
tions of sources in IRAC and MIPS images are less pre- 
cise due to the poorer angular resolution and larger pixel 
sizes in these bands. We selected the best astrometric 
matches to each Spitzer source from the objects identi- 
fied by SExtractor within a specified search radius, 9. 

4 http: / /ssc. spitzer. caltech.edu/dataanalysistools/tools/contributcd/ 
irac/epe/ 

5 http: / / ssc.spitzer.caltech.edu/dataanalysistools / tools /mopex / 

6 http: / / ssc.spitzer.caltech.edu / dataanalysistools / cookbook/ 
23/#.Toc256425880 

7 http:/ /ssc. spitzer. caltech.edu/irac/iracinstrumenthandbook/ 
IRAC_Instrumcnt_Handbook.pdf 
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To determine the best value of 9, we scrambled the RA 
of SExtractor sources and determined how many Spitzer 
sources were matched to a scrambled galaxy as a func- 
tion of 9. We found the best balance between purity and 
completeness for 9 ~ 1'.'25. This search radius yielded 
spurious matches for less than 2% of objects. The ac- 
tual contamination of our catalog will be much lower, 
because a Spitzer object with a spurious match will usu- 
ally be better matched to its "true" counterpart, which 
has a median match distance d = 0'.'4. The images used 
to perform the matching do not suffer from substantial 
confusion, even in the cluster centers, so erroneous pho- 
tometry due to overlapping sources is unlikely to present 
a problem. Further details of the visible image reduction 
were described by M06. 

2.3. Photometric Corrections 

We estimated the Galactic reddening toward each of 
the 8 clusters in ou r sample from the dust map of 
iSchlegel et al.l (|1998| ) and calculat ed extinction c orrec - 
tions assuming Ry = 3.1 and the iCardelli et al. (1989) 
redden ing law. The resolution of the ISchlegel et al.l 
(1998) dust map requires us to use a common extinction 
correction for all cluster members. However, Galactic cir- 
rus is apparent in some of our images, so this assumption 
is not always appropriate. This leads to additional un- 
certainty associated with the extinction corrections, but 
the total (visual) extinction toward our clusters is typ- 
ically less than 0.1 mags. The associated uncertainties 
are therefore small. For the clusters with the highest ex- 
tinctions (Abell 2104 and 2163, with A v = 0.73 and 1.1, 
respectively), variations in extinction across the cluster 
represent an important source of systematic uncertainty. 
We account for this by adopting a 10% uncertainty in all 
extinction corrections and propagating this uncertainty 
to the corrected magnitudes. In Abell 2163, for exam- 
ple, this corresponds to an uncertainty of 0.11 mags in 
the dereddened 1^-band magnitude. 

The raw fluxes measured from the MIR mosaics must 
be corrected for various instrumental effects, including 
aperture, array-location and color corrections, as de- 
scribed in the IRAC and MIP33 Instrument Handbooks. 
Aperture corrections are, in principle, required for all 
observations. In practice, even our smallest apertures 
(~ 7") are large enough that aperture corrections to 
visible- wavelength point sources are negligible. For MIR 
point sources, this is not the case. We apply aperture cor- 
rections from the IRAC Instrument Handbook appropri- 
ate for our redshift-dependent photometric apertures to 
the IRAC photometry. These corrections are not strictly 
appropriate due to the extended nature of our sources; 
however, we have chosen apertures that are large com- 
pared to the sources (~ 3x larger than the FWHM of 
the largest galaxies, see Section [Q]) . We therefore apply 
aperture corrections appropriate for point sources. 

We determined aperture corrections appropriate for 
our MIPS images by averaging a theoreticalpoint-source 
response function (PRF) from STinyTinfl with three 
bright, isolated point sources in the Abell 3125 and Abell 

8 http:/ /ssc. spitzer. caltcch.edu/mips/mipsinstrumenthandbook/ 
MIPS_Instrument_Handbook.pdf 

9 http:/ /ssc. spitzer. caltech.edu/dataanalysistools/tools/ 
contributed / general / stinytim / 



2104 mosaics. The PRFs of sources from the different 
clusters agree with one another and with the theoretical 
PRF to within a few percent over the range of aperture 
sizes relevant for our MIPS photometry. The dispersion 
between the individual PRFs at fixed aperture size pro- 
vides an estimate of the uncertainty on the corrections 
and is included in the 24/xto error budget. The MIPS 
images of the other clusters lack bright, isolated points 
sources with which to make a similar measurement, so 
we assume that the PRF appropriate for Abell 3125 and 
Abell 2104 gives reasonable aperture corrections for all 
clusters. This introduces some systematic error in our 
derived 24/zm fluxes, but the agreement of the observed 
PRFs of point-sources identified in Abell 3125 and Abell 
2104 with the theoretical PRF indicates that this uncer- 
tainty is small. 

The flatfield corrections applied to IRAC images by 
the automated image reduction pipeline are based on 
observations of the zodiacal background light, which is 
uniform on the scale of the IRAC field of view. It is 
also extremely red compared to any normal astrophys- 
ical source. The combination of scattered light due to 
the extended nature of the source and the color of the 
source illuminating the detector for the flatfield images 
results in different gains for point-sources and extended 
sources. It also requires an effective bandpass correction 
that varies with position on the detector. These effects 
can be corrected by applying a standard array-location 
correction image to a single IRAC image. For a mosaic, 
the magnitude of the required correction is significantly 
reduced by adding dithered images with different cor- 
rections at a given position on the sky. However, the 
residual effect can be a few percent or more depending 
on the number of overlapping IRAC pointings. We con- 
struct an array-location correction mosaic by co-adding 
the correction image for a single IRAC pointing shifted 
to the positions of each dithered image in the science mo- 
saic. We measure the required array- location corrections 
in the same apertures used to measure the IRAC fluxes. 

The Spitzer image reduction pipeline assumes a flat 
power-law SED to convert electrons to incident fluxes. 
Astrophysical sources typically do not show flat SEDs 
and therefore require color corrections to determine the 
true flux at the effective wavelength of a given band. This 
is especially important in star-forming galaxies, which 
show strong polycyclic aromatic h ydrocarbon (PAH ) 
emission features at 6.2 and 7.7/xm (Smi th et al.l 120071) . 
We determine color corrections to the measured fluxes 
from model SEDs (Section l3.1[) . We compute preliminary 
model SEDs for each cluster member from the photome- 
try with all other corrections applied. We then integrate 
the model SED across the various MIR bandpasses and 
determine the appropriate color corrections following the 
procedures outlined in the instrument handbooks. The 
color correction, K, applied to an IRAC source is given 

by, 

S{F v /F„ ){ V / Vo )-iRAv m 
S{v/v Q )-*R v dv 1 ' 

where F v is the model spectrum and R v is the re- 
sponse function of the detector in the appropriate chan- 
nel. The formalism for MIPS color corrections is similar 
but slightly more complicated; we refer interested read- 
ers to Section 3.7.4 of the MIPS Instrument Handbook. 
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Optical and MIR photometry for each cluster member 
after all relevant corrections have been applied are listed 
in Table 2. 

3. METHODS 

We wish to identify cluster members hosting AGNs, 
determine the AGN luminosities, examine the properties 
of AGN host galaxies, and determine whether they differ 
in any appreciable way from "normal" cluster galaxies or 
from their counterparts in the field. This requires that we 
distinguish cluster members from foreground and back- 
ground galaxies, fit model SEDs to the member photom- 
etry, and measure the rest-frame properties of the AGN 
host galaxies. We describe the model SEDs in Section 
13.11 Using these models, we calculate K-corrections to 
the measured fluxes, estimate stellar masses and SFRs 
for cluster member galaxies, an d identify AGN s. 

We use redshifts reported in iMartini et al.l (|2007l ) or 
extracted from the NASA Extragalactic DatabascFI to 
identify members of the galaxy clusters in our sample. 
We define a galaxy to be a cluster member if it falls 
within a circular field with radius, 



R < i?2oo = 1.7/T 1 Mpc 
[(I + z^ + Oa]- 1 / 2 



1000 km s- 1 



(3) 



wher e a is the cluster's velocity dispersion (jTreu et al.l 
[2001 . We also require that members have spectro- 
scopic redshift s within the ±3 ( 7 reds hift limits prescribed 
in Table 1 of IMartini et al.l (|2007l ). which were estab- 
lished using the bi weight velocity dispersion estimator of 
iBeers et al.1 (|1990T ). This criterion yields a sample of 1165 
cluster member galaxies. We eliminate many of these 
galaxies from our sample due to either limited photomet- 
ric coverage or, in a few instances, because the spectro- 
scopic redshifts in the literature are clearly in disagree- 
ment with the photometric redshifts obtained from the 
SED fits (Section [3T|). The final sample of "good" clus- 
ter members, those galaxies with detections in at least 
5 bands and with apparently reliable spectroscopic red- 
shifts, contains 488 galaxies. 

3.1. Model SEDs 

Assef et al. (2010; hereafter A10) constructed empiri- 
cal SED templates that can be used to determine photo- 
metric redshifts and K-corrections for galaxies and AGNs 
over a wide range of redshifts. The A10 templates include 
three galaxy templates (elliptical, spiral, and starburst or 
irregular) and a single AGN template, which can be sub- 
jected to variable intrinsic reddening. These templates 
were derived empirically across a long wavelength base- 
line (0.03-30/im), using 14448 apparently "pure" galax- 
ies and 5347 objects showing AGN signatures. We fit two 
independent model SEDs to the photometry of each clus- 
ter member using the published codes of A10. The first 
model included only the three galaxy templates, while 
the second also included an AGN component. The x 2 
differences between the two fits can be used to identify 
AGNs (Section Model SEDs for the M06 X-ray 

point sources included in our sample of "good" galaxies 

10 http:/ /nedwww. ipac.caltech.edu/ 



are shown in Figure [T] AGNs identified from their SED 
fits, but which have no X-ray counterparts, are shown in 
Figure[2] The fits to the X-ray point sources are represen- 
tative of the fit quality returned for all cluster members, 
while the fits to photometrically-identified AGNs are, on 
average, poorer. 

The model SEDs fit to 25 of the 488 spectroscopically- 
identified cluster members are poorly matched to the 
measured photometry (x 2 > 25). We determine photo- 
metric redshifts for all of the identified cluster members, 
and in cases where the measured photometric redshifts 
are more than 3tr away from the cluster redshift, we re- 
place the spectroscopic redshifts with photometric red- 
shifts and repeat the fit. In 11 cases, this procedure re- 
sults in substantial improvements to the fits (Ax 2 > 12, 
Xphoto-z < 4). This suggests that some galaxies in the 
sample have erroneous spectroscopic redshifts. One such 
object is an X-ray source, identified as AC 114-5 by M06. 
The redshift for this object was reported by Couch et al. 
(2001; their galaxy #365). The spectra used by these 
authors covered a relatively narrow wavelength range 
(8350 A < A < 8750 A) and had moderately poor S/N. 
We suspect that this combination of factors, in concert 
with a strong prior in favor of cluster membership in 
the presence of a putative Ha emis sion line at the cor- 
rect redshift, led iCouch et al.1 ((2001) to mis-identify the 
[On]AA4354 and [Om]AA4363 emission lines of a back- 
ground quasar at z — 0.988 as the [Nn]AA6548 and Ha 
emission lines, respectively, at the cluster redshift. Four 
of the 5 objects flagged as having erroneous re dshifts in 
AC 114 have redshifts from ICouch et al.l (|2001f ). Two of 
the four have redshifts from only one emission line, and 
we have confirmed that both objects with redshifts from 
multiple emission lines have plausible pairs of lines near 
the photometric redshifts. Furthermore, all of the ob- 
jects with apparently erroneous redshifts are quite faint, 
having V < 22, which makes acquiring high-S/N spectra 
difficult. Our identification of objects with discrepant 
photometric and spectroscopic redshifts as interlopers 
appears to be reliable, and we eliminate the associated 
galaxies from further consideration. The absence of AC 
114-5 from the X-ray AGN sample has important reper- 
cussions, which we discuss in Section |4j 

3.2. AGN Identification 

We consider AGNs selected based on their X-ray lumi- 
nosities, the shapes of their SEDs, or both. X-ray sources 
with Lx > 10 42 erg s _1 are unambiguously AGNs, but 
a number of processes can produce X-ray luminosities 
in the 10 40 - 10 42 erg s _1 range. These include LMXBs, 
HMXBs and a galaxy's extended, diffuse halo gas. The 
integrated X-ray luminosities of LMXBs and hot halo 
both correlate strongly with stell ar mass, as measured by 
the galaxy's -RT- band luminosity ()Kim fc Fabbianol 12004 
iSun et al.ll2007t) . a nd the luminosity f rom HMXBs cor- 
relates with SFR (jGrimm et al.ll2003D . These correla- 
tions allow us to predict the X-ray luminosity of a nor- 
mal galaxy using only parameters that can be measured 
fro m the model SEDs. Simila r analyses we r e perf ormed 
by iSivakoff et all (f2008h and I Arnold et all (|2009f), who 
used if-band luminosities measured from 2MASS pho- 
tometry rather than luminosities estimated from model 
SEDs. 
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Fig. 1. — Model SEDs for galaxies hosting M06 X-ray point sources. Bands shown are, in order of wavelength, B, V, R, I, [3.6], [4.5], 
[5.8], [8.0] and [24.0]. The panels are labeled with the names assigned by M06 in their Table 4. Objects also identified as AGNs from their 
SED fitting are labeled "IR." The heavy lines show the total model SED, while the solid, dotted, dashed and dot- d ashed lines show the A10 
AGN, elliptical, spiral and irregular templates, respectively. Not all components appear in all panels. See Section [XT] for further details. 



We measure if-band magnitudes from the model SEDs 
and determine SFRs from the K-corrected 8/im and 
24/xm luminosities of X-ray sources in each cluster. We 
use Lk and SFR in Eqns. 21 [5] and [6] to predict the ex- 
pected X-ray luminosities from t he host galaxies of X-ray 
point sources ident i fied by M06 (iKim &: Fabbianol 120041 : 
iGrimm et all 1200$ ISun et al.ll2007l respectively). The 
predictions for X-ray emission from a given galaxy due 
to LMXBs, HMXBs and the thermal halo are good to 
within ^0.3 dex and are given by, 



L x (thermal; 0.5-2 keV) 



2.5 x 10 39 erg s 



10 U L, 



1.63±0.13 



(6) 



L X (LMXB: 0.3 - 8 keV) = 
[(0.20 ± 0.08) x 10 30 erg s _1 ] 



J K,Q 



L X {HMXB) 



2.6 x 10 39 erg s" 1 



SFR 



M Q yr~ 



1.7 



(4) 



(5) 



where Lk and Lk s are the galaxy's luminosities in the K- 
and X s -filters. Each relation is given in slightly different 
energy ranges, none of which coincide with the range used 
by M 06. This probl e m is e specially severe for Eqn. [5j be- 
cause IGrimm et al.l ()2003[ ) take their X-ray fluxes from 
various sources in the literature without converting them 
to a common energy range. They claim that the resulting 
uncertainty is small because the scatter in the relation 
is much larger than the bandpass corrections. Fortu- 
nately, even if this were not the case, the HMXB contri- 
bution to the total predicted X-ray luminosities is small 
for the SFRs typical of cluster galaxies (< 10 M yr _1 ). 
The contribution from thermal emission to the soft X- 
ray luminosity can be significant, dominating the LMXB 
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Fig. 1. — Continued. The poor fit to AC114#5 indicates a bad spectroscopic redshift. 
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Fig. 2. — Model SEDs for objects identified as IR AGNs which are not also identified as X-ray AGNs. Line types and bandpasses shown 
are the same as in Figure Q] The object names indicated on each panel correspond to those in Table [T] See Section 13. II for further details. 
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component for L so f t > 6 x 10 erg s . This transi- 
tion luminosit y depends on the sp e cific f orm adopted 
in Eqn. H IMulchaev fc Jeltemal (|2010( ) found that 
Lx{corona) oc L^ 9±0A for fiel d galaxies , which differs 
sig nificantly from the results of Sun et all ()2007l ) . While 
the IMulchaev fc Jeltemal (|2010 ) relation is not strictly 
applicable to our sample, the difference between cluster 
and field galaxies suggests that the thermal X-ray emis- 
sion from a galaxy's halo depends on its environment. 
Such a variation introduces a systematic uncertainty in 
Lx(corona) of up to 0.8 dex at Lk = 4x 10 11 L q . Here- 
after we neglect this uncertainty, as its effect in a given 
cluster is impossible to quantify given the data presently 
available. 

We convert Eqns. 0][6]to determine luminosities in the 
soft X-ray (0.5-2 keV) and hard X-ray (2-8 keV) bands, 
assuming a T = 1.7 power law for the LMX B and HMXB 
relatio ns. We further assume that the iGrimm et al.l 
(2003) relation corresponds to luminosities in the 2- 
10 keV range and that the thermal emission from the 
kT = 0.7 keV halo gas is negligible in the hard X-ray 
band. The X-ray luminosities reported by M06 and our 
estimates of the systematic uncertainties in these lumi- 
nosities associated with the choice of energy correction 
factor (ECF) are shown in Figure [31 along with the pre- 
dicted luminosities from the host galaxies. Many of the 
reported point sources require an AGN component, but 
several of the M06 point sources have very massive host 
galaxies, and their observed fluxes may arise entirely 
from non-AGN sources. 

M06 selected 40 X-ray point sources with reliable de- 
tections above the extended emission from the surround- 
ing ICM (N count > 5). Of these 40 sources, they iden- 
tify 35 as probable AGNs. We have sufficient photome- 
try to construct reliable model SEDs for 35 M06 X-ray 
point sources. The remaining 5 M06 point sources either 




10 39 10 40 10 41 10 4Z 10 39 1 40 1 41 10 42 
Host X-ray Luminosity (erg s _1 ) 

Fig. 3. — Comparison of the X-ray luminosities of X-ray point 
sources from M06 (y-axis) to the predicted X-ray luminosities of 
their host galaxies (x-axis). Points show the measured luminosities, 
and the "tails" connect each source to the luminosity estimated by 
separating its 0.5 — 8.0 keV X-ray luminosity into soft and hard 
components using a T = 1.7 power-law. The length of the tail 
indicates how well the measured photon energies are described by 
a r = 1.7 power law, and consequently describes the systematic 
uncertainty on the quoted Lx. Long tails belong to objects poorly 
described by a T = 1.7 power law. Heavy lines mark the line of 
equality (Lx = L^ost), and the dashed lines show the ±0.7 dex 
scatter about the empirical relations used to pre dict the X-ray 
luminosity of a given host galaxy. See Section 13.21 for the method 
used to predict X-ray luminosities of normal galaxies. 

lack enough data to produce a reliable model SED or 
fall outside the i?-band field of view. We find that 23 of 
these 35 sources have X-ray luminosities more than la 
greater than the predicted host luminosity. Henceforth, 
we will call these objects X-ray AGNs. The systematic 
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flux error estimates in Figure [3] indicate that many X-ray 
AGNs have photon energy distributions that are poorly 
matched to the T = 1.7 power-law assumed by M06. 
Three such AGNs are close to the boundary separating 
probable AGNs from more ambiguous cases and have too 
large a soft X-ray flux compared to their hard X-ray flux 
to be consistent with a T = 1.7 power law. M06 did 
not correct for X-ray absorption, and in the cases where 
the ratio of soft to hard X-ray photons is too low for a 
r = 1.7 power law, absorption may explain the apparent 
discrepancy. However, objects whose soft X-ray fluxes 
are unexpectedly large compared to the total cannot re- 
sult from absorption. 

Many narrow- line Seyf ert 1 galaxies (NLS 1) show ex- 
cess soft X-ray emission ([Arnaud et al|[l985l ). However, 
only one X-ray source identified by M06 is a NLS1 (their 
Abell 644 #1), so the soft X-ray excess common to 
NLSls cannot explain the presence of excess soft X-ray 
emission in 13 X-ray sources with AGN-like luminosi- 
ties. Alternative explanations include soft X-rays arising 
from gas that is p hotoionized by an obscured AGN (e.g. 
IGhosh et al.ll2 0071. poor signal-to-noise in the X-ray, and 
thermal emission from hot gas. The ECF used to convert 
soft X-ray photons to incident fluxes for kT = 0-7 k eV 
thermal bremsstrahlung (assumed by iSun "eTaI1[200l is 
larger than the ECF for a T = 1.7 power law by ap- 
proximately 10%. This implies that two of the three sus- 
pect X-ray AGNs have luminosities sufficiently close to 
the threshold that they may reasonably be mis-classified 
galaxies. This yields a possible contamination in the X- 
ray AGN sample of approximately 10%, which is com- 
parable to the estimated contamination of the IR AGN 
sample (see below). 

In comparison to our sample of 23 X-ray AGNs from 
a parent sample of 35 X-ray point sources with com- 
plete photometry, M06 found that 35 of their 40 point 
sources had X-ray luminosities consistent with AGNs. 
The larger fraction of AGNs reported by M06 may be 
attributed to their use of Lx~Lb relations, which show 
larger scatter than the K-band relations. We also intro- 
duce some uncertainty by estimating Lk from the model 
SEDs, but this uncertainty is small (~ 10%) compared to 
the scatter in the L^-Lk relation. An additional differ- 
ence is that M06 considered the two luminosity compo- 
nents separately and did not compare their sum to the 
measured lum i nositi es, T his was done s ubsequ ently by 
iSivakoff et all (|2008t ) and lArnold et a\\ (|2009l ) in their 
studies of AGNs in low-redshift groups and clusters of 
galaxies. Their analyses are much closer to our method, 
and their samples included some of the clusters in our 
sample (Abell 3128, 3125 and 644). 

An alternative method to identify AGNs is to use 
the distinctive shape of their SEDs, particularly in 
the MIR (e.e iMarconi et all [200l iStern et all [20051 : 
iRichards et all [20061 : A101 This approach can identify 
AGNs behind gas column densities large enough to ob- 
scure even the X-rays emitted by an AGN. Such an AGN 
sample has very different selection criteria and biases 
than an X-ray selected sample, and combining the two 
results in more complete AGN identification. 

We identify AGNs from their SEDs by comparing the 
goodness-of-fit of two sets of model templates. The first 
set uses only the normal galaxy templates. The other 
also includes the AGN template. We determine whether 




Fig. 4. — Likelihood ratio (p) distributions for fits to model galax- 
ies with no AGN component after photometric errors have been 
added. Each panel shows the distributions resulting from exam- 
ining 10,000 normal galaxies. Galaxies with p = 1 are well-fit by 
the three normal galaxy templates and do not require an AGN 
component. An object with p = would be perfectly fit by the 
4-template model SED and have X^ a i+AGN = dashed line 

indicates th e se lection threshold, p m ax used to identify IR AGNs. 
See Section [32] for further details. 

a given galaxy requires an AGN component in its model 
SED by applying a threshold on the likelihood ratio, p, 



P : 



exp[-x 2 (gaQ/2] 
cxp[- X 2 {gal + AGN)/2] 



(7) 



where ^(gal) and ^{gal + AGN) are goodnesses-of-fit 
for a model with only the A10 galaxy templates and for 
a model that includes an additional AGN component, 
respectively. AGNs are those objects whose p is smaller 
than a pre-determined selection limit, p m ax, established 
by Monte Carlo simulations of normal galaxies. 

We created artificial galaxy photometry to determine 
an appropriate p m ax by combining the three galaxy tem- 
plates of A10 in proportions that reflect the template 
luminosity distributions in real cluster members. We in- 
troduced Gaussian photometric errors comparable to the 
photometric uncertainties in our real data (0.07 mag) to 
the fluxes given by the model SEDs. We also allowed 
occasional catastrophic errors of up to 0.3 dex. The 
artificial galaxy photometry did not include upper lim- 
its, which we also neglected when constructing model 
SEDs of real galaxies. We fit the artificial galaxies with 
two models. The first model excluded the AGN compo- 
nent from the fit, while the second component included 
it. The likelihood ratio distributions computed from the 
goodness-of-fit results for the two different models are 
shown in Figure 2] These distributions show the proba- 
bility that a pure galaxy will be erroneously classified as 
an AGN due to the presence of photometric errors. The 
similarity of the different distributions, even based on 
only 4 photometric bands, indicates that a single p m ax 
can be used to select AGNs from among all galaxies in 
our sample. 

We also identify AGNs based on the F-statistics of the 
two model SED fits described above. Figure [5] shows 
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the F-statistic as a function of ^(gal) for X-ray AGNs 
selected using Figure [3j AGNs selected using likelihood 
ratios, and "normal" cluster members. The F-statistic is 
given by, 



F 



Ax 2 /2 



Xl(gal + AGN) 



(8) 



where A% 2 is the (absolute) change in the total % 2 after 
introducing the AGN component to the fit. In addition 
to the galaxies that are well-fit by the galaxy-only model 
and not substantially improved by the addition of an 
AGN component, there are objects with large x 2 (gaZ) 
but small F, and objects with large F but small X 2 (ff a 0- 
Neither of the latter two categories contain objects likely 
to be AGNs from the point-of-view of the model SEDs. 
The most luminous X-ray AGNs have both large F and 
large x 2 (ff°0- These are clearly identified as AGNs by 
the model SEDs, and less luminous X-ray AGNs can be 
found with increasing density toward the normal galaxy 
locus at the origin of Figure [5] The dotted and dashed 
lines in the Figure correspond to the p < p ma x selec- 
tion boundaries for N=6 and N=9 flux measurements, 
respectively. Some objects above the N=9 line are not 
selected as IR AGNs because they fail a cut on the over- 
all goodness-of-fit, which requires X^{gal + AGN) < 5. 
We could define an AGN selection region in Figure [5l but 
due to the non-uniformity of our photometric data, this 
would result in different effective cuts in A\ 2 between 
different clusters and between objects in individual clus- 
ters. Furthermore, we find that only 3 AGNs identified 
using likelihood ratios fall into the suspect part of Figure 
4 with F « 1. This level of contamination (~ 10%) is 
consistent with the estimated purity of the X-ray AGNs, 
which we deem to be acceptable. Therefore, for the rest 
of this work, we rely on the more simplistic likelihood 
ratio threshold to identify AGNs. 

Likelihood ratio selection of AGNs using SED fitting is 
most sensitive to the shape of the MIR SED, so we refer 
to AGNs so identified as IR AGNs. We find 29 IR AGNs 
using a selection boundary at the 99.8% confidence in- 
terval of the merged p distribution (p max = 1.5 x 10~ 3 ). 
Table [I] lists both X-ray and IR AGNs, their luminosi- 
ties, and the basic parameters of their host galaxies. IR 
AGN selection re covers 5 of 7 AGN s (71%) identified via 
the Stern wedge (jStern et al.l 120051 see Figure H]) and 8 
of the 23 X-ray AGNs. The galaxies in the Stern wedge 
that are not selected from their SED fits fall just in- 
side the boundary of the wedge, so they may be normal 
galaxies shifted into the wedge by photometric errors. 
IGoriian et all (120081) find that 35% of X -ray sources m 
the Bootes field of the NO AO Deep Wide Field Survey 
(fx > 8 x 10~ 15 erg s _1 cm~ 2 ) with detections in all 4 
IRAC bands fall outside the Stern wedge, and Figure 14 
of A10 shows that a substantial fraction of the point- 
source (luminous) AGNs in their sample fall outside the 
wedge as well. Given the high luminosities in both of 
these samples, it is perhaps not surprising that most of 
the lower-luminosity AGNs common in galaxy clusters 
fall outside the Stern wedge. For p max = 1.5 x 10 ~ 3 and 
the size of our sample (488 galaxies), we expect on aver- 
age one false-positive AGN identification and 3 or fewer 
false-positives at 98% confidence, implying > 90% purity 
in our IR AGN sample. 




Fig. 5. — Application of F-statistics for AGN identification based 
on fits of galaxy-only and galaxy-plus- AGN model SEDs to mea- 
sured photometry. Red triangles show AGNs identified using the 
p threshold shown in Figure [4] (IR AGNs). Open blue squares 
show X-ray AGNs, and solid black squares show "normal" cluster 
members. The dotted and dashed curves show the p thresholds for 
objects having N=6 and N=9 flux measurements, respectively. Ob- 
jects above their corresponding selection boundaries are identified 
as AGNs, provided that they pass a \i cu t- 

We estimate the completeness of the IR AGN sample 
as a function of the reddening of the AGN template and 
luminosity using Monte Carlo simulations. We construct 
model AGN SEDs by injecting an AGN component with 
some luminosity and reddening into artificial galaxy pho- 
tometry, which we generate using the Monte Carlo tech- 
niques described above. We estimate the completeness 
of the IR AGN sample from the fraction of such AGNs 
recovered. The completeness depends strongly on the lu- 
minosity of the AGN component. We only reliably iden- 
tify AGNs with L bot > 7 x 10 10 i Q . The completeness 
depends only weakly on E(B — V). There are measur- 
able differences only for AGNs with E(B — V) > 2. For 
our observed wavelengths, AGN identification depends 
most strongly on the shape of the MIR SED, which is 
insensitive to modest amounts of reddening. The full de- 
pendence of completeness on Lf, i and E(B — V) is listed 
in Table H 

We caution that both our AGN identification and the 
analysis below were conducted using the fixed AGN tem- 
plate derived by A10. While this template is dominated 
by luminous AGNs, AGNs of all luminosities were used 
in its construction, and in some sense it represents the 
optimal median AGN SED. There is some evidence that 
AGNs with low Eddington ratios (Lboll^Edd) are sys- 
tematically weaker in the UV and the MIR than higher 
Lboi/LEdd AGNs. This appea rs to become important at 
L bo i/L E dd « 1CT 3 (|Hol 120081 ). However, the UV weak- 
ness of such objects remains a subject of debate (e. 
lHol[T999l 120081 IDudik et aLlUM lEracleous et"aUl201C 
and the SEDs of AGNs appear to all be quite similar out 
to A 20pm, even in AGNs w ith accretion rates as low 
as L ho i/L Edd w 10" 3 (|Hol 120081 Figure 7). Furthermore, 
the variable reddening of the AGN component allowed 
by the models can account for differing UV/visible flux 
ratios, making the AGN component of the model SEDs 
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Fig. 6. — Positions of both IR AGNs (red triangles, upper panel) and X-ray AGNs (blue pentagons, lower panel) on the Stern ct al. (2005) 
AGN selection diagram. The dashed, green trapezoid marks the Stern wedge. Filled symbols mark AGNs whose colors were determined 
using only measured magnitudes; open symbols show colors determined using model SEDs to estimate magnitudes in missing IRAC bands. 
No model colors were constructed for normal galaxies. Passive galaxies are located in the lower-left corner, normal star-forming galaxies 
are found in the center, and (U)LIRGs occupy the upper-right part of the galaxy sequence. The error bars in the upper-left show the color 
uncertainty for a typical cluster member, including AGNs. 



flexible enough to mimic AGNs with a wide variety of 
Eddington ratios. 

Intrinsic variations in the AGN SED are one possi- 
ble cause of the absence of an important AGN compo- 
nent in the SEDs of many X-ray AGNs, despite their 
similar distributions in Lj, ; (Section [4]). Another possi- 
ble explanation is that the nuclear MIR emission from 
many X-ray AGNs is overwhelmed by star-formation 
in their host galaxies. We find that X-ray AGNs with 
Lx > 10 42 erg s" 1 that are also identified as IR AGNs 
have no measurable star-formation, while those not iden- 
tified in the IR have (SFR) = 0.3 M Q yr" 1 . This may 
be a selection effect, since nuclear MIR emission is not 
subtracted before computing SFRs in galaxies not iden- 
tified as IR AGNs. However, it appears that the balance 
between SFR and nuclear emission is an important fac- 
tor in determining whether a given X-ray source will be 



identified as an IR AGN. 

Also of concern is the MIR emission exhibited by 
some normal galaxi es which is clearly not associated with 
star-f ormation (e.g. iVerlev et al.ll2009t Kelson fc Holdenl 
120101 ). The strength of the diffuse interstellar dust 
emission relative to star formation varies from galaxy 
to galaxy depending on the populatio ns of AGB stars, 
which can produce and heat dust (jKelson fc Holdenl 
|2010|), an d field B-stars (including HB stars), which pro- 
duce UV light that can both heat du st grains and excit e 
PAH emission in the diffuse ISM fe.g. lLi fc Drainell2002l) . 
These effects could mimic the presence of an AGN, par- 
ticularly in passively-evolving galaxies, which the A10 
templates predict should decline strictly as a vF v oc j/ 2 5 
power-law. Given the limited data available to constrain 
MIR emission not associated with either an AGN or a 
star-forming region and the as-yet uncertain magnitude 
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of the associated variations, we neglect any potential 
effects on our AGN identification. However, potential 
sources of MIR emission not accounted for by the A10 
templates, especially emission from dust heated by old 
stars in passive galaxies, remain a potentially important 
systematic uncertainty. 

3.3. Stellar Masses 

Stellar population synthesis modeling provides a means 
to es timate stellar masses in the absence of detailed spec- 
tra. iBell fc de Jond (|2001|) construct model spectra of 
galaxies for a wide variety of stellar masses, SFRs, metal- 
licities and stellar initial mass functions (IMFs) to con- 
vert colors to mass-to-light ratios (M/L). Their models 
assume a mass-dependent formation epoch with bursty 
star-formation histories, which is appr opriate for the 
spiral galaxies they study. Figure 9 of IBell &: de Jond 
(2001) makes it clear, however, that their results also ro- 
bustly estimate M/L for passively evolving galaxies. In 
fact, the scatter about the mean M/L tends to decrease 
for redder systems because the stochasticity of the star- 
formation history becomes less important in galaxies that 
experienced their last burst of star-formation in the dis- 
tant past. 

Bell & de Jong provide a table of coefficients (a\,b\) 
relating M/L for a galaxy to its color, 



log 10 (M/L\) = a x + b x x color 



(9) 



where color is measured in the bands for which a\ 
and 6a were determined. We adopt the coefficients 
appropriate for Solar m etallicity computed with the 
iBruzual fc Chariot! (|2003l ) population sy nthesis code and 
the sc aled Salpeter IMF suggested by Bell & d e Jongl 
(|2001D . who report that a modified Salpeter IMF with 
total mass M' — 0.7Ms a ipeter yields the best agreement 
with the Tully-Fisher relation. Once we select an ap- 
propriate (a\,b\) pair, it is straightforward to compute 
stellar masses from the visible photometry. However, we 
must first subtract the AGN component of the model 
SED in sources identified as IR AGNs before computing 
colors. 

The uncertainty introduced by the AGN subtraction is 
a combination of the fractional uncertainty in the contri- 
bution of the AGN template to the model SED, which is 
determined by the fit, and the uncertainty in the AGN 
template itself. To measure the uncertainty in the tem- 
plate, we examined 1644 luminous quasars with spec- 
troscopic redsh ifts from the AGN and G alaxy Evolution 
Survey (AGES; lKochanek et alJlin prenl ) and determined 
the variation in their measured photometry about their 
best-fit model SEDs. Using these measurements, we con- 
structed an RMS SED for AGNs and averaged it across 
each of the bandpasses we employ. The uncertainty in the 
AGN correction resulting from intrinsic variation about 
the AGN template is 10% except at 24/im, where there 
are too few z — quasars to make a meaningful compar- 
ison. The uncertainty in the AGN correction at 24/xm is 
therefore large, but it can be constrained by the relatively 
good agreement of the 8/xm and 24/ito SFRs (Figure [7]). 

In galaxies with no genuine nuclear activity the AGN 
template can correct for variations in stellar populations 
relative to the templates, intrinsic extinction, or errors in 
the measured photometry. Subtraction of the AGN com- 
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Fig. 7. — Comp arison of S FRs determined using the IRAC and 
MIPS relations of lZhu et all pOOl l. The solid lines in both panels 
denote equality. The bottom panel shows the fractional residuals 
between the two star-formation indicators. The MIPS SFRs are 
biased high with respect to IRAC SFRs by 0.1 dex, which is sig- 
nificant at 2.4<t but smaller than the intrinsic scatter (0.2 dex). 

ponent under these circumstances would result in under- 
estimated stellar masses and SFRs, while failure to sub- 
tract the AGN component in a genuine, low-luminosity 
AGN would cause the measured SFRs of their host galax- 
ies to be biased toward higher values. However, the am- 
biguity between a genuine, low-luminosity AGN and an 
apparent AGN component introduced to correct for pho- 
tometric errors (Section [32]) renders any attempt to sub- 
tract the AGN component in such cases suspect. There- 
fore, in normal galaxies and in X-ray AGNs not identified 
as IR AGNs, no AGN correction is applied. We accept 
the inherent bias to avoid introducing ambiguous AGN 
corrections, which would be much more difficult to inter- 
pret. 

The IBell fc de Jond (|2001[ ) calibrations are reported 
for rest-frame colors, so we need K-corrections for each 
cluster member to convert the measured magnitudes to 
the rest frame. We calculate the K-corrections from the 
model SEDs returned by the A10 fitting routines. Uncer- 
tainties on K-corrections cannot be directly determined 
from the uncertainties in the model components because 
K-corrections depend non-linearly on these uncertainties. 
Therefore, we recombine the components of each model 
SED in proportion to the uncertainties in their contri- 
butions to the total model flux. This results in a series 
of temporary model SEDs. We then calculate the K- 
corrections implied by these temporary model SEDs and 
measure their dispersions to estimate the uncertainties in 
the K-corrections returned by the original model SED. 

The systematic uncertainty on the stellar masses cal- 
culated from Eqn. [9] can be estimated by comparing 
the fiducial masses with masses derived using different 
assumptions. We estimate the typical systematic uncer- 
tainty in stellar mass, listed in Table [TJ to be 0.2 dex. 
These uncertainties are derived by measuring the differ- 
ence between the fiducial masses and those determined 
using coefficients appropriate for the Pe gase popula- 
tion synthesis models with a Salpeter IMF. lConrov et all 
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(2009) studied the ability of different models to repro- 
duce the observed colors of stellar populations in glob- 
ular clusters and found that systematic uncertainties on 
stellar masses derived from population synthesis codes 
typically reach or exceed 0.3 dex. 

3.4. Star- Formation Rates 

We measure SFRs from our AGN-corre cted MIR 
photo metry using the empirical relations of iZhu et aLl 
(2008), which have been determined for both the IRAC 
and the MIPS 24^m bands using the same calibra- 
tion sample. While the contribution of the stellar con- 
tinuum to the observed 24/um luminosity is negligible, 
the Ray leigh- Jeans tail of the stellar continuum emission 
can make an important contribution to the integrated 
flux at 8/im, especially in galaxies with the low SFRs 
typical in clusters. The method used to subtract this 
contribution is an im portant sy s tema tic uncertainty in 
the SFR calculation. IZhu et ail (|2008f ) assume that the 
contribution of the stellar continuum at 8/j,m can be de- 
scribed by L s v tellar ( 8um) = 0.23 2 ^(3. 5^777), as derived 
from the models of IHelou et al.l (|2004l ). Under this as- 
sumption, IZhu et all (|2008l) derive luminosity-SFR rela- 
tions appropriate for a Salpeter IMF, 



SFR(M Q yr- 1 ) 



SFR(M Q yr" 1 ) = 



isLf, ust (8fim) 
1.58 x 10 9 L o 

vL v (2Ap,m) 
7.15 x 10 8 L Q 



(10) 



(11) 



where L^ us *(8/im) is determined by subtracting 
F s J ellar (8/im) from the the measured 8fim luminos- 
ity. Simoes-Lopes et al. (in preparation) find that 
L s u tellar {8nm) = 0.269X^(3.5^771) provides a better 
estimate for their sample of nearby, early-type galaxies 
with no dust and conclude tha t the d ifference in their 
result compared to IHelou et al.l (|2004l ) is due to metal- 
licity. Another important systematic uncertainty in 
SFRs derived from PAHs i s the dependence of the PAH 
abundance on metallicity (Calzetti et al. 2007J) because 
lower metallicity systems have fewer PAHs and therefore 
weaker 8fim emission at fixed SFR. This second effect is 
negligible for the high-mass — and therefore metal-rich — 
galaxies we consider. We neglect both metallicity- and 
mass-dependent effe cts for the r emain der of our analysis. 
Instead, we follow IZhu et all (|2008[ ) and assume that 
L s v tellar (8pm) = 0.232L !/ (3.5/jm). We derive SFRs from 
Eons. [TU1 and [TT1 For galaxies having measurable (> 3cr) 
SFRs from both IRAC and MIPS, we take a geometric 
mean of the two; otherwise, we use whichever SFR 
measurement is available. The resulting SFRs for AGNs 
are summarized in Table [TJ 

Eq uat ions ITU1 and [TT1 were derived using the extinction- 
corrected Ha luminosity of the associated galaxies. The 
MIPS SFR determined from Eqn. [TT] for a galaxy with 
vL v = 7.15 x 1 9 L^ is re 0.6 dex lar ger than the SFR de- 
rived from the Calzetti et al.l (|2007l ) relati on, which was 
calibr ated using the Paa emission line. ICalzetti et al.l 
(|2007f) used the Starburst99 IMF, and after accounting 
for this difference, the resulting discrepancy is reduced 
to 0.4 dex. The choice of SFR calibration therefore rep- 
resents an important systematic uncertainty in the mea- 
sured SFRs. The total systematic uncertainty in SFR 



is indicated by the significant scatter (0.2 dex) and the 
small but marginally significant offset (0.1 dex) between 
the IRAC and MIPS SFRs in Figure[7] Since the offset is 
smaller than both the scatter about the line of equality 
and the systematic u ncertainty when comparing to the 
ICalzetti et al.l (|2007f ) result, we neglect it below. How- 
ever, we caution that there remains a ~ 15% uncertainty 
in our results associated with the discrepancy between 
the IRAC and MIPS SFR indicators. 

4. RESULTS 

We identify 29 IR AGNs with likelihood ratios p < 
Pmax- We also confirm the presence of AGNs in 23 X- 
ray point sources whose X-ray luminosities significantly 
exceed the luminosities expected from their host galax- 
ies. Surprisingly, the X-ray and IR AGN samples are 
largely disjoint: only 8 AGNs appear in both. Only the 
more luminous IR AGNs appear in the X-ray AGN sam- 
ple and vice-versa. While it is not surprising for faint 
X-ray AGNs to drop out of the IR AGN sample, the ab- 
sence of X-ray emission associated with many IR AGNs, 
which require a moderately luminous AGN for a reliable 
detection, is unexpected. This may indicate either dif- 
ferent selection biases in the two methods or genuine, 
physical differences between the AGNs selected by these 
techniques. 

4.1. Bolometric AGN Luminosities 

In order to conduct a meaningful comparison of X-ray 
and IR AGNs, we need to place them on a common lumi- 
nosity system. The most obvious choice is the bolometric 
AGN luminosity (£b z), which also allows us to examine 
black hole growth rates. 

The A10 AGN template provides a natural means 
of determining the bolometric luminosity (Ltoi) for IR 
AGNs, but the MIR luminosity in the template comes 
from reprocessed dust emission, which would result in 
double-counting the U V emission from th e disk for 
AGN s viewed face-on ([Marconi et al.l [2004!. hereafter 
M04; iRichards et al.l I2006T) . We instead determine Lboi 
using a piecewise combination of the AGN model SED 
and three power-laws. We integrate the unreddened A10 
AGN template from Lya to 1/_ito, shortward of which 
the template becomes uncertain due to absorption by the 
Lya forest. We estimate the X-ray luminosity by inte- 
grating a r = 1.7 power law from 1-10 keV. We estimate 
the extreme ultraviolet (EUV) luminosity by integrating 
L v cx v- a °* from A = 1216A to 1 keV. The slope of 
the EUV SED (a ox ) is given by Eqn. 2 of IVignali et al.l 
(2003]), 



= 0.1 log 



L„(2500A) 



erg s 



- 1.32 



(12) 



with L iy (2500A) taken from the AGN template SED. Fi- 
nally, we eliminate reprocessed emission from dust by 
assuming F v cx v~ 2 for 1/im < A < 30/im, following 
M04. 

To correct the X-ray luminosities of X-ray AGNs to 
bolometric luminosities, we fit a power-law to the mea- 
sured I/ X (0.3 - 8 keV) and L hoi of the 8 IR AGNs identi- 
fied separately in X-rays. A least-squares fit to the total 
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FlG. 8. — Comparison of our estimated relation between Lf, i and 
Lx (black solid line) with that found by M04 (red dashed line). The 
measured X-ray (0.3-8 keV) and bolometric luminosities of AGNs 
identified by both the X-ray and IR selection criteria are shown by 
black points. L bo i is derived by integrating the A10 AGN model 
SED component from 1216A to 1/im and assuming a declining con- 
tinuum with F u oc v~' 2 for A > 1 /im. Lx is determined assuming 
a r = 1.7 power law. See Section 14.21 for further information. 
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X-ray and AGN luminosities yields, 

log[L x (0.3 - 8 keV)] = (0.9 ± 0.2) log 

+ (41.4 ±0.2) (13) 

where Lboi is the bolometric AGN luminosity integrated 
from 10 keV to 30/Ltm. The AGNs used to determine 
Eqn. [13] show a scatter of 0.4 dex about the best-fit 
relation (Figure H]). Figure [8] suggests that the slope 
returned by the fit may be strongly influenced by the 
highest-luminosity AGN. However, a fit to the other 7 
AGNs returns an identical slope (0.9±0.5), so Eqn. [TBI is 
not significantly biased by the highest-luminosity object. 
The luminosity dependence of the bolometric corrections 
(BCs) derived from the fit is therefore robust. The slope 
is also consistent, within large statistical uncertainties, 
with the luminosity-dependence derived by M04. 

The BCs derived from Eqn. [T3] are fairly crude. For ex- 
ample, the fit does not account for uncertainties on Lx 
or Lb i. It also ignores upper limits, which will lead it 
to over- predict the true Lx at fixed Lboi- M04, by con- 
trast, provide luminosity-dependent BCs in several en- 
ergy ranges that account for X-ray non-detections (their 
Eqn. 21). We convert their BCs to 0.3-8 keV assum- 
ing r = 1.7 and estimate the expected X-ray flux from 
our IR AGNs. The predicted X-ray fluxes exceed those 
estimated using Eqn. I13[ which we know over-estimates 
the intrinsic Lx~Lboi relation, by 0.7 dex or more. This 
might result if the M04 SED is a poor match to the A10 
AGN template. M04 deter mine their X-ray BCs using 
the a ox relation derived by iVignali et al.l ([2003) for a 
sample of SDSS quasars, including broad-absorption line 
quasars (BALQSOs). Given that our Lboi calculation is 
insensitive to the absorption in BALQSOs, it is possible 
that the M04 BCs over-estimate L x at fixed L„(2500A) 
when applied to our sample. In order to produce consis- 



tent results for the X-ray and IR AGNs, we therefore use 
the BCs implied by Eqn. [H rather than the M04 BCs, 
despite the large uncertainties associated with Eqn. 1131 

4.2. X-ray Sensitivity 

With the Lx-Lboi relation provided by Eqn. [T3l we 
can determine whether the X-ray non-detection of many 
IR AGNs results from some intrinsic difference between 
the two classes of AGNs or if it is merely a result of the 
sensitivity of the X-ray images used by M06. Eqn. [13] 
predicts that 9 (5) IR AGNs with no X-ray detections 
should be more than a factor of 3 (5) brighter than the 
faintest point source in their parent clusters (M06). The 
M04 BCs produce more X-ray flux at fixed bolometric lu- 
minosity than Eqn. [T3]and yield 13 (12) IR AGNs with 
significant X-ray non-detections with the same flux lim- 
its. The lack of detectable X-rays from many IR AGNs 
is consequently easier to explain if we use Eqn. [13] rather 
than the M04 relations to predict their X-ray luminosi- 
ties. 

The minimum detected flux in a given cluster may not 
always be a fair representation of the sensitivity for a 
given IR AGN due to variations in the Chandra effec- 
tive area with off-axis angle. However, the magnitudes 
by which many IR AGNs in AC 114 exceed the mini- 
mum detected flux, sometimes more than a factor of 5, 
suggest that these AGNs should have been detected if 
they obeyed the Lx~Lboi relation of Eqn. [13] The non- 
detection of many IR AGN s in X-rays is quali tatively 
consistent with the results of Hickox et al.l (|2009| ) , whose 
IR AGN selection relied upon the Stern wedge, and who 
found many strong IR AGNs that could not be identi- 
fied in X-rays. At least some of the "missing" IR AGNs 
could be highly obscured. An intervening absorber with 
Nh = 10 22 cm -2 would reduce the observed 0.5-2 keV 
flux by a factor 3, which is sufficient to explain many of 
the missing IR AGNs. The missing AGNs could also re- 
sult from the large scatter about the mean Lv(2500A)- 
a ox relation. The AGN with the most significant X- 
ray non-detection exceeds the minimum reported flux by 
a factor of 7 , whic h can be explained by Ao OI ~ 0.4. 
IVignali et al.l ([2003) report a large intrinsic scatter about 
their best- fit relation, and the combination of this scatter 
with in situ absorption could mask moderately luminous 
AGNs from detection in X-rays. 

Finally, at least one IR AGN (A1689 #109) appears 
to be absent from the M06 sample due to X-ray vari- 
ability rather than as a result of absorption, intrinsic X- 
ray faintness, or shallow Chandra imaging. This object 
is moderately luminous {Lboi = 2.1 x 10 10 L Q ), AGN- 
dominated (Jagn = 0.95), falls firmly in the middle of 
the Stern wedge, and is very robustly detected by our 
likelihood ratio selection (p = 4 x 10~ 77 ). Nevertheless, 
there is no X-ray point source associated with this object 
in the Chandra image employed by M06. In a more re- 
cent observation (Chandra Obs ID 6930, PI G. Garmire), 
A1689 #109 is associated with an X-ray point source 
far brighter than the X-ray sources reported by M06. 
It therefore seems likely that the IR AGNs that require 
the most extreme values of a ox could be accounted for 
by variability rather than by systematically weak X-ray 
emission compared to their visible- wavelength luminosi- 
ties. 
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Fig. 9. — Cumulative stellar mass, SFR and sSFR distributions of 
the X-ray and IR AGN samples compared to the distributions for 
all cluster members. Neither of the AGN samples show any signif- 
icant differences in either Mr or SFR compared to the full sample 
of cluster members, nor does the merged AGN sample. However, 
IR AGN hosts have higher sSFR than both X-ray AGN hosts and 
normal galaxies at 99% confidence, despite their similarities in M* 
and SFR. 

4.3. Host Galaxies 

We determine stellar masses and SFRs for AGN host 
galaxies after subtracting the AGN component from the 
SED. This introduces some additional uncertainty in the 
resulting masses and SFRs beyond the original photo- 
metric uncertainties, as discussed in Section 13.31 The 
uncertainty in the AGN contribution to the measured 
MIR fluxes can prevent detection of low-level star for- 
mation in IR AGNs. The SFR distribution among IR 
AGNs is therefore biased toward high SFR. 

Figure [9] shows the results of comparing galaxies host- 
ing different types of AGNs to one another and also to 
cluster galaxies as a whole. The stellar mass and SFR 
distributions of galaxies hosting X-ray and IR AGNs 
show no measurable differences with the distributions of 
all cluster members. Merging the X-ray and IR AGN 
samples likewise yields no measurable difference. How- 
ever, the hosts of IR AGNs have high specific SFRs 
(sSFR) compared to the hosts of X-ray AGNs and to 
all cluster members at 98% and 97% confidence, respec- 
tively. The difference between the sSFRs of X-ray AGN 
hosts and the full galaxy sample is not significant. How- 
ever, X-ray AGN hosts appear to have lower sSFRs than 
the average galaxy in Figure |9l whi ch is consistent wit h 
previous results using field galaxies (Hickox et al. 2009). 
We must also consider the effect of non-detections on the 
measured distributions. Many of the IR AGN hosts have 
upper limits on SFR that are smaller than the SFRs of 
the X-ray AGN host galaxies with the lowest measurable 
SFRs. Therefore, if the IR AGN hosts had a distribution 
of SFRs similar to the X-ray AGN hosts with measurable 
star-formation, star-formation would have been detected 
in most IR AGN hosts. This indicates that uncertain- 
ties in the AGN corrections alone cannot account for the 
higher sSFRs among IR AGN hosts. 

The IRAC color-color diagram (e.g. IStern eFaIll2005h 



probes of the nature of AGN host galaxies independent 
of their model SEDs by identifying the dom inant source 
of their MIR emission (jDonlev et all 120081 ). The MIR 
colors of X-ray and IR AGNs before their AGN com- 
ponents are subtracted are compared to all cluster mem- 
bers in Figure |6] Galaxies hosting AGNs have unremark- 
able [5.8] — [8.0] colors but do not extend as far to the 
red as normal galaxies, which indicates that AGNs are 
seld om found in starbu rsts or luminous infrared galax- 
ies (jDonlev et all 12008?) . AGN hosts also show redder 
[3.6] — [4.5] colors than typical for a red sequence galaxy, 
which may indicate a contribution of hot dust to the 
4.5/im continuum. The colors of AGN hosts, especially 
IR AGN hosts, are influenced by the AGN continuum, 
but tests using the AGN and spiral galaxy templates in- 
dicate that only galaxies in the Stern wedge have more 
than 50% of their IRAC fluxes contributed by the AGN 
component. A two-dimensional KS-test confirms that, 
after excluding objects in the Stern wedge, the IRAC 
colors of both X-ray and IR AGNs differ from normal 
galaxies at > 99.9% confidence, and the absence of X- 
ray AGNs among the most vigorously star- forming galax- 
ies (those with the reddest [5.8] — [8.0] colors) is consis- 
tent with earlier indications that X-ray AGNs avoid the 
blue cloud in visible co lor-magnitude diag rams (CMDs; 
iSchawinski et all l2009t IHickox et ail 120091 ) . The distri- 
bution of X-ray AGNs in Figure [6] also appears to be 
consistent with the results of IGoriian et al.l (|2008D . who 
found that 16.8 ±0.3% of X-ray-identified AGNs outside 
the Stern wedge had very red [5.8] — [8.0] colors consis- 
tent with vigorous, on-going star-formation. We found 
this population to be 20 ± 6% among our X-ray AGNs. 

The visible CMD provides a means to estimate the 
nature of galaxies in the absence of measurable star- 
formation at MIR wavelengths. Figure [10] shows the 
CMD for each cluster after the AGN component has been 
subtracted. The fraction of cluster members hosting an 
X-ray AGN peaks on the red sequence, and the probabil- 
ity that the X-ray AGN hosts are drawn from the parent 
cluster population is less than 10 -3 . This contrasts with 
AGN hosts in the field, where the X-ray AGN fraction 
typically peaks in the green valley (IHickox et alJ 12009b 
ISchawinski et alj|2009b iSilverman et al.ll2009l henceforth 
S09) for AGNs identified using either X-ray luminos- 
ity or emission-line diagnostics. IR AGN hosts, both 
in our sample of clust er AGNs and in the field sample 
of IHickox et al.1 (|2009h . conspicuously avoid the red se- 
quence. Like the difference between X-ray AGN hosts 
and the parent cluster population, this result is signif- 
icant at > 99.9% confidence. This indicates that the 
IR AGN sample has at most limited contamination by 
MI R-excess early-type galaxies of the sort studied by, 
e.g. iBrand et all ([2009). Galaxies hosting IR AGNs in 
clusters also show an important difference compared to 
their counterparts in the field. While only 1.5% of field 
galaxi es hosting the IR AGNs studied by IHickox et al.l 
(2009) had 01 (u — g) colors redder than the median of 
the red sequence, more than 20% of IR AGNs in clusters 
have visible colors redder than the red sequence in their 
parent clusters. 

We examined the SDSS g — r colors of very red galaxies 
(XV — R)rest- frame > 0.8) in Abell 1689, which has the 
largest number of such objects, and found that most also 
appear red in SDSS colors. The most notable exception 
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Fig. 10. — Optical color-magnitude diagrams showing the 
spectroscopically-confirmed member galaxies (black squares), X- 
ray AGNs (blue hexagons) and IR AGNs (red triangles) of each 
cluster. The contribution of the AGN component to the model 
SED has been subtracted from the IR AGNs, leaving estimated 
host-galaxy colors and luminosities. Typical uncertainties on colors 
and absolute magnitudes are approximately 0.1 mag. See Section 
I4.3l for further discussion. 



is Abell 1689 #192, which we have identified as an IR 
AGN, which suggests that its colors may change due to 
AGN variability. The qualitative agreement between the 
colors of very red galaxies in Figure [TOl and their g — r col- 
ors from SDSS suggests that these objects are genuinely 
unusual and not the result of photometric errors. These 
galaxies also show substantial reddening of the AGN tem- 
plate in their A10 fit results, with (E(B-V)) = 0.4 and a 
trend for higher E(B — V) in galaxies with redder colors 
at 97% confidence. These results suggest that the un- 
usually red galaxies in Figure QJj] experience significant 
internal extinction that is not present in most galaxies. 

Since the AGN component of the SED fit may account 
not only for a true AGN contribution but also for intrin- 
sic variations about the normal galaxy templates, some 
or all of these very red AGNs, which represent approx- 
imately 1/3 of our IR AGN sample, may not be true 
AGNs. However, fewer than half (7/17) of objects with 
(V — R)rest- frame > 0.8 arc identified as IR AGNs; this 
implies that IR AGNs must differ from normal galaxies 
not only in the visible but also in the MIR, and MIR 
fluxes are practically immune to extinction. Therefore, 
most of the IR AGNs identified in this region of color- 
magnitude space cannot be false-positives selected due 
to their unusual visible colors but must have genuine nu- 
clear activity contributing to their SEDs. 

4.4. Accretion Rates 

We use the bolometric luminosities of both X-ray and 
IR AGNs to measure the growth of their black holes and 
compare the black hole growth to the assembly of stellar 
mass in their host galaxies. The accretion rate of a black 
hole can be generically written as, 
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Fig. 11. — Comparison of accretion rates (Mbh', left panel) and 
SMBH growth rates relative to their host galaxies (Mbh / SFR; 
right panel). The various distributions compare different meth- 
ods of estimating Mbh- directly from the SEDs of IR AGNs (red 
dashed), applying ad hoc BCs derived from AGNs with both X- 
ray and IR identifications (solid black), and applying the M04 BCs 
(blue dotted). M04 BCs return significantly lower Mbh than the 
other two methods, which are consistent with one another. Short 
black and red arrows mark upper limits for X-ray and IR AGNs, 
respectively. The dashed vertical line mar ks the ratio required t o 
maintain th e z = MBH-Mij U i ge relation (Marconi & Hunt 2003). 
See Section l4.4l for more on the various BCs. 

where L^oi is the bolometric luminosity and e is the 
efficiency of conversion between the rest mass energy 
(Mc 2 ) of the accreted material and the energy radiated 
by the black hole. We assume e = 0.1, appropriate for 
a thin accretion di sk around an SMBH with moderate 
spin (jThornd Il974) and determine L^oi as described in 
Section O 

The accretion rates derived from Eqn. [JJ] for the X- 
ray and IR AGN samples are shown in Figure [TT] The 
left panel suggests that X-ray and IR AGNs have similar 
accretion rates, and a KS test reveals that there is no sig- 
nificant difference between the two samples. This is sur- 
prising, since we would naively expect that the difference 
between X-ray and IR AGNs might be due to different 
dependence of X-ray and IR selection techniques on lumi- 
nosity. Instead, the right panel of Figure [TTJ shows that 
the X-ray and IR AGN samples have (M B h/SFR) = 
3 x 10~ 3 and (Mbh /SFR) = 2 x 10~ 3 , respectively. 
These ratios are comparable to the mean Mbh /Mbuige 
in the local universe (2 x 10~ 3 . IMarconi fc Huntl [2003h . 
which indicates that the SMBHs in cluster AGNs are ac- 
creting at approximately the rate required to maintain 
the z = MsH-Mbuige relation. However, this is likely 
an artifact of our SFR detection thresholds, as the accre- 
tion rates of these objects are not large enough to pro- 
duce outliers on the Msn-Mbuige relation in a Hubble 
time. 

Figure[T^compares black hole accretion rates with host 
mass and sSFR. We find no significant correlation be- 
tween Mbh and sSFR, nor do we find a correlation of 
Mbh with stellar mass among the X-ray AGN sample. 
However, Mbh correlates with stellar mass among IR 



18 



Atlee et al. 



0.1 



c-0.01 



10- 




10" 12 10" 11 10" 10 10 
sSFR (yr" 1 ) 



10 1 



M, 



10 1 
(Me) 



Fig. 12. — Relationships of black hole accretion rates (Mbh) to 
stellar masses and sSFRs. Black points and arrows s how Mgg 
inferred from X-ray luminosities using BCs from Eqn. 1131 while 
red points and arrows mark Mgu inferred from integrating model 
SEDs. Stellar masses and SFRs are measured from SEDs and in- 
clude the entire galaxy, not just the spheroidal component. 

AGNs at 99.5% confidence, weaking to 98% confidence 
among the merged AGN sample. This correlation may 
be related to the ability of more massive cluster members 
to retain more cold gas. 

Figure [13] shows the relationship between black hole 
growth and stellar mass assembly in AGN host galax- 
ies. The correlation of Mbh with SFR is extraordinarily 
strong (> 99.9% confidence), and both X-ray and IR 
AGNs appear to follow the same relation, with SFR cx 



M 



0.46±0.06 
BH 



. iNet zcr (2009) studied emission line selected 
AGNs from SDSS and also found a tight correlation be- 
tween SFR and AGN luminosity across nearly 5 dex in 
L bo i. However, their SFR-M BH relation (SF R oc M%%) 
is stee per than ours at 5.7a. Furthermore, iLutz et al.l 
(2010) performed a stacking analysis of X-ray identified 
AGNs at z ~ 1 and found no measurable correlation of 
SFR with L boi for AGNs with L 2 -iokcV < 10 44 erg s -1 . 
However, t he millimete r- brigh t, optically luminous QSOs 
studied by ILutz et al.l ( 20081 ) appear to be con sistent 
with both iNetzerl (|20091 ) and ILutz et al.l (|201JJ.jrhe 
qualit ative similarity of our results to those of INetzerl 
(2009) and ILutz et ahl (|2010| ) suggests that we are see- 
ing the same underlying relationship. That both X-ray, 
IR and emission-line selected AGNs appear to show the 
same general trend toward higher Mbh in hosts with 
higher SFR suggests that accretion rates in all of these 
objects are set by the size of the global cold gas reser- 
voir. Such a relationship is also predicted theoretically 
as a result of large-scale dynamical instabilities, which 
drive col d gas to the centers of galaxies where it can be 
accre ted (jKawakatu & Wada 2008; Hopk ins fc Quataertl 
12010( 1. However, the quantitative discrepancies between 
the various observational signatures of star-formation 
and gas accretion indicate that further work on the rela- 
tionship between these phenomena is needed. 

Figure [13] also compares star-formation and black hole 
growth among our AGN sample with the median ratio 



found by S09 and the ratio needed to maintain the z = 
MsH-Mbuige relation. In some cases Mbh I 'SFR falls 
more than a dex below the ratio reported by S09 for 
field galaxies at z ~ 0.8 and more than 0.3 dex below 
the rate needed to maintain the local MBH~M bu i ge rela- 
tion. However, if we consider AGN hosts with no measur- 
able star- formation, the disagreement in Mbh/M* be- 
tween the cluster AGNs we measure and the field AGNs 
of S09 becomes far less pronounced. The upper limits 
in Figure [13] fill in much of the empty space between 
the S09 median relation and the cluster AGNs with 
measurable star-formation, but the fraction of galax- 
ies with Mbh / SFR < 2 x 10~ 3 is larger in Figure 
Hp than in Figure 13 of S09 (7/39 versus 9/67). This 
difference grows (7/27) if we consider only AGNs with 

Mbh < 10~ 2 Mq yr -1 , which is below the luminos- 
ity limit of the S09 sample. However, even the differ- 
ence between the low-luminosity subsample and the S09 
result is not statistica lly significant (90% confidence). 
iSilverman et al.1 (|2009f ) project the evolution in the me- 
dian SFR of their AGN sample to z = and find that 
it agrees with the SFRs measured in Type 2 AGNs with 
log(L [0///] ) > 40.5 in SDSS. The median z = 0.2 SFR 
for the S09 AGN hosts is SFR « 0.5 M yr" 1 , which is 
comparable to our detection threshold. As a result, the 
AGNs measured in Figure [13] are more comparable to a 
high-SFR subsample of the S09 AGNs. However, there is 
no significant difference in the Mbh / SFR of high-SFR 
versus low-SFR AGNs in S09. We therefore concluded 
that the ratio of Mbh to SFR our sample of low-z clus- 
ter AGNs is consistent with the ratios observed in high-z 
AGNs in the field. 

4.5. Radial Distributions 



iMartini et all (|2007f) found that luminous (L x > 
10 42 erg s _1 ) X-ray AGNs were more centrally concen- 
trated in R/R200 than normal cluster members at 97% 
confidence. After pruning the AGN sample of suspect 
rcdshifts and applying improved K-corrections, we as- 
semble the radial distributions of our AGN samples in 
FigureHU FiguresfMb and [14b. which consider the X-ray 
and IR AGN samples, respectively, have slightly different 
distributions of parent galaxies. This is because Spitzer 
pointings cover only the fields around X-ray sources iden- 
tified by M06 and not the full Chandra field of view. The 
IR AGNs are selected from the cluster member catalog 
after SED fitting has been performed, so the radial dis- 
tribution of IR AGNs is guaranteed to be unbiased with 
respect to the cluster galaxy sample we used above, while 
X-ray AGNs must be compared to the distribution of all 
galaxies within the Chandra footprint. These different 
selection footprints lead to the different radial distribu- 
tions shown in the solid red and black lines in Figure 
114b . The difference is not significant, however, and has 
no impact on our conclusions. 

We have determined that the host galaxy of the X-ray 
point source identified as the cluster AGN AC114-5 by 
M06 had an erroneous spectroscopic redshift reported in 
the literature (see Section l3~T1 and Figured]). Our SED 
fitting indicates that this source is a background QSO at 
Zphot ~ 0.99. Without this object, which is located at a 
projected distance R/R200 ~ 0.2 from the center of AC 
114, the significance of the difference between the lumi- 
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Fig. 13. — Relationship between black hole growth and star-formation in our AGN sample compared to the AGN sample from COSMOS 
examined by S09. Green circles and arrows mark the S09 AGNs. All other symbols are the same as in Figure | 12 | Lines mark th e Mbh /SFR 
relation measured by S09 (solid) and the ratio required to produce the z = MBH~Mbulge relation (Marconi & Hunt 2003; dashed). Our 
SFRs and those reported by Silverman ct al. (2009) are galaxy-wide SFRs rather than bulge SFRs. 



nous X-ray AGN and control samples drops to 89% confi- 
dence with a luminosity-selected control sample and 92% 
confidence with a mas s-selected cont r ol sam ple. Consis- 
tent with the results of lMartini et al.1 (|2007t ) , we also find 
no significant difference between the radial distribution 
of the full X-ray AGN sample compared to the distribu- 
tion of clust er members as a who le. 
Following iMartini et al.l (|2009Q . we also try a redshift- 



dependent luminosity threshold yM^ cut = Mjj(O) + 1 — 

in place of a fixed value. The galaxy and AGN 

samples selected using this criterion show no significant 
differe nces in their R/R200 distributions. IMartini et al.1 
(2009) chose this evolving threshold to select a sample of 
passively-evolving galaxies at fixed stellar mass. A mass 
threshold (M* > 3 x 1O 1O M0) appropriate for an ellipti- 



cal galaxy at z = with Mr = M^ cut again yields no 
measurable difference between the radial distributions of 
X-ray AGNs and all cluster members. We conclude that 
the radial distributions of both X-ray and IR AGNs in 
galaxy clusters are consistent with the distribution of 
cluster members, although the agreement between clus- 
ter members and IR AGNs is much better than between 
cluster members and X-ray AGNs. 

5. DISCUSSION 

Identifying AGN from their X-ray emission is widely 
considered to be among the most robust means of se- 
lecting AGNs (e.g. lUeda et~aLll2003l S09, A10), because 
the measured hard X-ray luminosity of a given AGN is 
largely insensitive to absorption if Nh < 10 cm~ 2 . Fur- 
thermore, the fraction of Compton-thick AGNs (Njj > 
10 24 cm~ 2 ) is small, with 10% or less of all cosmic 
black hole growth taking place in Compton-thick systems 
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R / R 200 

Fig. 14.— Distributions of X-ray AGNs, IR AGNs and all 
cluster galaxies within their parent clusters. Panel a: luminous 
(Lx > 10 42 erg s -1 ; dot-dashed) and all X-ray AGNs (dashed) 
compared to all cluster members with M* > M C ut = 3 X 10 10 Mq 
(solid). Panel b: luminous (Lf, ; > 10 1 " Lq; dot-dashed red) 
and all IR AGNs (dashed red) compared to cluster members with 
Spitzer coverage and M* > M cu t (solid red) and the unbiased sam- 
ple of galaxies inside the Chandra footprint (solid black). 

(jTreister et all 120091 ). Alternatively, AGNs can also be 
robustly identified from their UV continuum emission af- 
ter it lias been absorbed by dust and re-emitted in the 
MIR. If these techniques are similarly immune to the ef- 
fects of absorption they should yield very similar AGN 
samples. Instead, we find that at most 15% of AGNs 
in galaxy clusters are identified by both X-ray and MIR 
techniques. 

Furthermore, it is clear that this dichotomy does not 
result solely from the relative luminosities of X-ray and 
IR AGNs. The IR AGN sample contains 5-9 objects 
that should have been detected in X-rays if their SEDs 
were similar to those AGNs identified using both selec- 
tion methods. The most prominent of these is Abell 1689 
#109, which has L^oi ~ 8 x 10 45 erg s -1 but was not de- 
tected in the Chandra image used by M06 to identify 
X-ray AGNs in Abell 1689. This AGN appears quite 
prominently in a subsequent Chandra image, indicating 
that its initial non-detection was most likely the result 
of X-ray variability. This example demonstrates that the 
absence of detectable X-ray emission from an AGN can- 
didate, even a fairly luminous one, does not necessarily 
preclude the presence of an AGN. However, Abell 1689 
#109 is not typical. The IR AGNs with significant X- 
ray non-detections are not necessarily the most luminous. 
Instead, they reside in the clusters with the deepest X- 
ray images. Indeed, all of the X-ray non-detections in 
AC 114 that fall within the Chandra image footprint are 
predicted to be at least 3 times brighter than the faintest 
reported X-ray point source. As a result, at least some of 
these non-detections could indicate contamination of the 
IR AGN sample by one or more of the effects discussed 
in Section l3~2l e.g. intrinsic variation in the AGN SED or 
dust heating by AGB carbon-stars. More observational 
and theoretical work on the dust emission in old stel- 
lar populations are required before the potential of these 



sources of MIR emission to mimic an AGN-like SED can 
be quantified. 

In the absence of detailed, calibrated models for "con- 
tamination" of MIR emission by old stars, we assume 
that this component is negligible. This implies that X- 
ray selection alone can miss a large fraction of moderate- 
to-low luminosity AGNs. This could have important im- 
plications for studies o f star-formation i n clu s ters using 
MIR luminosities (e.g. Sainton ge et alj [20081 : iBai et all 
I2009t IGeach et al.ll2003T This is especially important 
if authors assume that AGNs can always be identified 
with X-rays alone or that the MIR emission from galaxies 
with X-ray excesses is always dominated by AGN emis- 
sion. These assumptions imply that any MIR emission 
not associated with an X-ray AGN must be powered by 
star-formation and that no MIR emission from a galaxy 
hosting an X-ray AGN can be powered by star-formation. 
Our results indicate that these assumptions may lead au- 
thors to overestimate the number of cluster galaxies with 
vigorous star-formation and to underestimate the num- 
ber with moderate star-formation. Therefore, additional 
tests for AGN are needed to correctly interpret the MIR 
luminosities of cluster galaxies. 

A difference between X-ray- and MIR-selected AGN 
samples also appears among field samples, which consist 
of more luminous AGNs than the o nes we study and us e 
a different MIR selection method ([Hickox et all 12009). 
The color distributions of IR AGNs selected using differ- 
ent techniques also differ from one another, but it is clear 
that galaxies hosting AGNs identified from their X-ray 
emission are dissimilar from galaxies hosting AGNs iden- 
tified in the MIR. Most notably, IR AGN hosts have sig- 
nificantly higher sSFRs than the average cluster galaxy, 
while there is no significant difference between the sSFRs 
of X-ray AGNs and the cluster population as a whole. 
Since SFR correlates well with cold gas mass, higher sS- 
FRs among IR AGN host galaxies suggests these galaxies 
have a larger fraction of their baryons in cold gas than 
X-ray AGN hosts. However, the differences discussed in 
Section r4.3l are determined only for galaxies with measur- 
able star-formation. Several IR AGNs are found in host 
galaxies that have both visible and IRAC colors consis- 
tent with passively-evolving stellar systems. 

The tight correlations between accretion rates of both 
X-ray and IR AGNs with SFR in their host galaxies sug- 
gests that the two classes are fueled by the same mech- 
anism and are therefore fundamentally similar. Subject 
to the caveat described above, the larger sSFRs found 
in IR AGN hosts might explain the apparent dichotomy 
of the two AGN classes despite their physical similar- 
ity. Larger gas fractions in IR AGN hosts could lead to 
larger average column densities in IR AGNs, depressing 
Lx/Lboi in these systems. The presence of at least 5 of 
the 8 IR AGNs with X-ray counterparts on the red se- 
quence, where there is little cold gas to participate in 
X-ray absorption, tends to support this scenario (Figure 
[TO)) . If cold gas fractions of AGN host galaxies influence 
the detectability of X-ray AGNs, this might also explain 
the dearth of X-ray AGNs in the green valley in clusters 
compared to the field. The X-ray AGNs in our sample 
are weaker than the AGNs usually studied in field galaxy 
samples, and a modest cold gas reservoir in green valley 
galaxies could more easily absorb enough X-rays from an 
AGN with Lx = 10 41 erg s _1 to make it undetectable. 



Cluster AGNs 



21 



Doing the same for an AGN with Lx = 10 43 erg s 1 , 
whi ch is more typical f or the field samples studied by, 
e.g. iHickox et all (|2009l ) and S09, would require a larger 
gas column. 

Just over half (58%) of the M06 X-ray point sources 
have detectable hard X-ray emission, and therefore many 
AGNs near the Chandra detection limits could be hidden 
by a sufficiently large absorbing column. Only 3 of the 
9 IR AGNs in AC 114 whose bolometric luminosities im- 
ply that they should have been detected in X-rays, but 
were not, would remain detectable in the soft X-ray band 
behind a gas column with Nh = 10 22 cm~ 2 . This col- 
umn density is large for Type I AGNs, but i t is not un- 
usual for Type II AGNs observed in X-rays (jUeda et al.l 
2003). Furthermore, X-ray and IR AGNs seem to obey 
the same relationship between SFR and accretion rate 
in AGN hosts whose SFRs are measurable. This is con- 
sistent with the hypothesis that the apparent dichotomy 
between X-ray and IR AGNs is false, and the shape of 
an AGN's SED depends strongly on the amount of ab- 
sorbing material between us and the central black hole. 

The scenario we propose, in which absorption by cold 
gas in the host galaxy is responsible for the absence of de- 
tectable X-ray emission from IR AGNs, is consistent with 
the differences we find between the two samples. How- 
ever, verifying that absorption by the host ISM is indeed 
the cause of this observed difference will require deeper 
X-ray observations to detect X-ray counterparts and es- 
timate absorption columns. If this can be accomplished, 
the presence of spectral signatures of X-ray absorption 
would confirm that the host galaxy is responsible for hid- 
ing some IR AGNs from X-ray detection. 

6. CONCLUSIONS 

We have used Spitzer imaging of galaxy clusters to 
identify AGNs and to measure the masses and star- 
formation rates of their host galaxies. We find that AGNs 
identified by this technique have very little overlap with 
AGNs identified in X-rays. We compared the host galax- 
ies of AGNs identified using the two methods and deter- 
mined that, while their masses and SFRs are indistin- 
guishable, IR AGNs reside in galaxies with higher sSFRs 
than both X-ray AGN hosts and the parent sample of 
cluster galaxies. The hosts of X-ray AGNs have sSFRs 
that are somewhat lower than but consistent with the 
sSFRs seen in cluster galaxies as a whole. The difference 
between X-ray AGN hosts and normal cluster galaxies is 
significant only when comparing their positions in visi- 
ble color-magnitude and MIR color-color diagrams. X- 
ray AGN hosts are rarely found in the regions of both 
diagrams associated with vigorous star-formation. 

We also find that accretion rates of both X-ray and IR 
AGNs correlate strongly with SFR in their host galaxies. 
This suggests that X-ray and IR AGNs are physically 
similar and are fueled by the same mechanism. We hy- 
pothesize that the larger sSFRs seen in IR AGN hosts in- 
dicate larger cold gas fractions in these galaxies, and sug- 



gest that this could account for the apparent dichotomy 
between X-ray and IR AGNs. A moderately large cold 
gas column density of 10 23 cm -2 could suppress the X- 
ray emission from the IR AGNs enough that we would 
be unable to detect them. The presence of IR AGNs 
but not X-ray AGNs in galaxies with very red optical 
colors, indicative of strong absorption, lends credence to 
this hypothesis. It might also be verifiable directly by 
deep X-ray observations of either AC 114 or Abell 1689 
to search for X-ray emission from IR AGNs and to de- 
termine if such X-ray emission shows evidence for ab- 
sorption intrinsic to the host galaxy. For example, the 
most luminous IR AGN with no X-ray counterpart in 
Abell 1689 could be detected by Chandra with S/N = 3 
per resolution element at 4 keV — the energy cutoff for 
objects with Ah = 10 23 cm -2 — in 160 ks. This would 
allow a crude model spectrum to be constructed and the 
intrinsic absorption column to be measured. Finally, we 
have obtained NIR spectra of several IR AGN in Abell 
1689, which we will examine for high-ionization emission 
lines that would unambiguously indicate the presence of 
an AGN. 

Following iMartini et al.1 ([2007), we compared the ra- 
dial distributions of AGNs and all cluster members. We 
eliminated one AGN with a spectroscopic redshift from 
the literature that incorrectly identified a background 
quasar as a cluster member. Without this object, the 
significance of their result that luminous X-ray AGNs 
(Xx > 10 42 erg s _1 ) are more concentrated than clus- 
ter members as a whole is reduced to ~ 90% confidence. 
While this result is no longer significant, it would be 
worthwhile to extend the present sample using archival 
Chandra imaging of additional clusters to either confirm 
or refute that X-ray luminous AGNs are more concen- 
trated than the galaxy populations of their parent clus- 
ters. It is unlikely, however, that a similar exercise using 
IR AGNs would yield a positive result, as the radial dis- 
tribution of IR AGNs agrees very closely with the distri- 
bution of cluster galaxies. 
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TABLE 1 
Cluster Sample 



Cluster 


z 

(1) 


(km s- 1 ) 
(2) 


^members 

(3) 


IRAC AOR(s) 
(4) 


vL v (, s (8//m) Limit 
(10 42 erg s- 1 ) 
(5) 


MIPS AOR(s) 
(6) 


vL v ^ b s {2A[im) Limit 
(10 41 erg s- 1 ) 
(7) 


Abell 3128 


0.0595 


906 


83 


25410816 


0.54 


25411072 


2.6 


Abell 3125 


0.0616 


475 


25 


25409792 


0.58 


25410048 


2.6 


Abell 644 


0.0701 


952 


9 


25409280 


1.0 






Abell 2104 


0.1544 


1242 


74 


25411328 


1.2 


25411584 


1.9 


Abell 1689 


0.1867 


1400 


160 


4754176, 14696192, 14696448, 14696704, 


1.3 


4770048. 4769792, 


2.8 










14696960, 14697216, 14697472, 25411840 




19042304, 19042048 




Abell 2163 


0.2007 


1381 


27 


25412352 


1.8 


25412608 


3.4 


MS 1008.1-1224 


0.3068 


1127 


68 


25410304 


0.81 






AC 114 


0.3148 


1388 


159 


4756480, 12653824, 25412864 


1.0 


4773888, 4774144, 25413120 


2.2 



Note. — Summary of clusters included in the analysis and the observations contributing to the MIR mosaic images of each cluster. The extra line beneath Abell 1689 contains 
additional AORs that do not fit on a single line. (1) Redshifts from Martini et al. (2007) f determined using the biweight estimator of Beers et al. (1990). (2) Velocity dispersions of 
cluster members estimated by Martini et al. (2009) using the biweight measure of Beers et al. (1990). (3) Total number of galaxies with both MIR and .R-band coverage identified as 
cluster members by Martini et al. (2007) or extracted from the literature using their redshift limits. (4) Astronomical Observation Request (AOR.) numbers of Spitzer observations used 
to contruct IRAC mosaics. (5) The minimum detectable observer-frame Sfim luminosity in each cluster, derived from the 3cr lower limit on measurable flux in a "typical" part of the 
8/im mosaic image. Due to the variable coverage across the cluster, lower luminosites are detectable in some cluster members than in others. (6) AORs used to contruct the 24p,m 
mosaics. (7) 3c" lower limits on detectable 24/jm luminosities. These are derived in a similar manner to the IRAC limits in column (5) and have the same caveats. 



TABLE 2 
Cluster Member Photometry 



Name 


RA 


Dec 


B 


V 


R 


I 




F„(3.6/im) 


F„(4.5inn) (mjy) 




F„(8.0iim) (mjy) 


F„(2iiim) 


















(mjy) 


(mJy) 


(mJy) 


(mjy) 


(mjy) 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




(8) 


(9) 


(10) 


(11) 


(12) 


a3128-001 


03:30:37.7 


-52:32:57 


17.88 ± 0.08 


16.81 ±0.06 


16.21 ±0.06 






1.37 ±0.08 


0.90 ± 0.08 


0.66 ±0.17 


< 0.77 


< 0.84 


a3125-001 


03:27:20.2 


-53:28:34 


17.87 ±0.12 


16.89 ±0.09 


16.29 ±0.09 








1.37 ±0.14 




9.49 ± 0.85 


8.68 ± 2.28 


a644-011 


08:17:39.5 


-07:33:09 


17.29 ±0.12 


16.69 ±0.09 


16.20 ±0.09 






2.61 ±0.21 


2.55 ±0.24 


3.62 ±0.48 


5.15 ± 0.62 




a2104-001 


15:40:07.6 


-03:17:06 


20.59 ±0.16 


19.34 ±0.12 


18.62 ±0.12 


18.05 ± 


0.11 


0.18 ±0.02 


0.13 ±0.01 


< 0.08 


< 0.19 


< 0.07 


al689-004 


13:11:29.5 


-01:20:27 


17.87 ±0.14 


16.48 ±0.10 


15.75 ±0.10 


15.16 ± 


0.09 


2.89 ± 0.26 


2.11 ±0.19 


1.21 ±0.11 


0.87 ±0.08 


< 0.29 


a2163-001 


16:15:25.8 


-06:09:26 


20.70 ±0.19 


19.69 ±0.16 


18.88 ±0.14 






0.10 ±0.01 


0.08 ±0.01 


< 0.07 


< 0.21 


< 0.10 


msl008-001 


10:10:34.1 


-12:39:52 


22.58 ±0.12 


21.11 ±0.09 


20.10 ±0.08 


19.41 ± 


0.08 


0.08 ±0.01 


0.07 ±0.01 


0.04 ±0.01 


< 0.05 




acl 14-001 


22:58:52.3 


-34:46:47 




22.40 ±0.12 


21.73 ±0.10 


21.40 ± 


0.11 


0.00 ± 0.00 


0.00 ±0.00 


< 0.01 


< 0.02 





> 
o 

2? 



Note. — Visible and MIR photometry for a small selection of example galaxies. The full table is available from the online version of the journal. (1) The name of this object, constructed 
from a shorthand of its parent cluster and the order in which each object appears in the list of cluster members extracted from NED. (2-3) Positions of this object in J2000 coordinates, 
as derived from the i?-band images. (4-7) Visible photometry for each object, where detectable, in Vega magnitudes. Fluxes are measured in the i?-band Kron-like aperture. Objects with 
no quoted magnitudes in a given band have either no coverage or no detection in that band. No upper limits are quoted. (5-8) MIR fluxes measured in ff-band Kron-like aperture. Where 
appropriate, 3<r upper limits on measured MIR fluxes, derived from the appropriate uncertainty mosaic, are given. Galaxies with no quoted upper limit for a given band have no coverage in 
the corresponding image. 



to 

00 



24 



Atlee et al. 



TABLE 1 
Identified Active Galactic Nuclei 



Name 


Martini Name 


RA 


Dec 


Hoi (lO^erp s x ) 


L x (10 4i er<; s i ) 


M*(10 ll> Mq) 


SFR (M© yr- 1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


a3128-004 


a3 128-9 


03:30:39.3 


-52:32:05 


1.5 ±0.5 


3.8 


3.2 ± 0.9(stat) ± 1.9(syst) 


< 0.04 


a3128-012 


a3128-6 


03:30:17.3 


-52:34:08 


1.0 ±0.3 


22.4 


1.9 ± 0.6(stat) ± l.l(syst) 


< 0.35 


a3128-092 


a3128-2 


03:29:41.4 


-52:29:35 




22.4 


1.1 ± 0.2(stat) ± 0.6(syst) 


< 0.44 


a3125-044 


a3125-5 


03:27:05.0 


-53:21:41 




7.2 


7.5±1.6(stat)±4.5(syst) 


2.7 ±0.4 


a644-011 


a644-l 


08:17:39.5 


-07:33:09 


10.4 ±1.4 


28.2 


2.6 ± 0.9(stat) ± 1.6(syst) 


< 1.38 


a644-024 


a644-2 


08:17:48.1 


-07:37:31 




4.5 


7.6±1.8(stat)±4.5(syst) 


1.2 ±0.4 


a2104-024 


a2104-4 


15:40:14.0 


-03:17:03 




7.2 


2.0 ± 0.6(stat) ± 1.2(syst) 


< 0.50 


a2104-040 


a2104-6 


15:40:03.9 


-03:20:38 




7.2 


12.4 ± 3.8(stat) ± 7.4(syst) 


1.9 ±0.6 


a2104-046 


a2104-5 


15:40:19.5 


-03:18:24 




36.3 


1.7±0.5(stat) ± l.O(syst) 


< 0.83 


a2104-051 


a2104-2 


15:40:16.7 


-03:15:07 


4.1 ±0.8 


18.2 


< 2.3 


< 16300.00 



Note. — Brief sample table summarizing AGNs identified either by their X-ray luminosity or their SED shapes. The full table is av ailable from the 
electronic edition of the journal. (1) The name of this object in Table 77. (2) The name given to the X-ray source bv lMartini et all f2006|) . (3-4) Position 
of this AGN in J2000 coordinates, as derived from the i?-band image. (5) The bolometric luminosity derived by integrating the direct component of the 
AGN cont ribution to the mode l SED. These luminosities arc quoted only for IR AGNs. (6) Rcst-framc X-ray luminosities in the 0.3-8 keV band from 
Table 4 of [Martini et al. (2006). X-ray luminosities arc given only for X-ray AGNs. (7) Stellar mass derived using the M/L coefficients appro priate for 
a solar metallicity galaxy with a scaled Salpeter IMF and applying the Bruzual & Chariot population synthesis model (Bell Sz dc Jong 2063, Table 4). 
Systematic errors arc derived by applying the M/L coefficients appropriate for a Salpeter IMF and the Pegase population synthesis model. Upper limits 
are given at 3cr of the statistical error only. (8) SFR derived cither from the Sp,m luminosity, the 24fim luminosity or by taking the geometric mean of 
the two, depending on the measurements available. Uncertainties include only statistical errors, and upper limits arc quoted at 3er in the more sensitive 
of the 8fim and 24/^m bands. 



TABLE 2 
IR AGN Selection Efficiency 



Lagn : 


-100 -> -1.0 


-1.0 -> -0.4 


-0.4 -> 0.2 


.2 -> .8 


0.8 -> 1.5 


1.5 -> 100 


E(B -V) 














0-0.075 


6.9 ±0.3% 


8.1 ±0.3% 


17.8 ± 0.4% 


45.5 ± 0.9% 


72.9 ± 1.9% 


83.8 ± 4.6% 


0.075-0.15 


7.0 ± 0.4% 


8.3 ± 0.3% 


17.8 ± 0.4% 


45.6 ± 0.6% 


75.5 ± 2.0%% 


83.8 ± 5.0% 


0.15-0.3 


7.0 ± 0.3% 


8.3 ± 0.3% 


17.1 ±0.4% 


44.7 ± 0.8% 


74.2 ± 1.7% 


84.3 ± 4.2% 


0.3-0.4 


6.7 ± 0.5% 


7.5 ±0.4% 


18.2 ± 0.6% 


46.2 ± 1.3% 


77.2 ± 2.7% 


86.8 ± 6.9% 


0.4-0.6 


7.2 ± 0.4% 


8.4 ±0.3 


17.8 ± 0.5% 


49.2 ± 1.0% 


77.0 ±2.1% 


85.8 ± 5.1% 


0.6-1.0 


7.2 ± 0.4% 


8.4 ± 0.3% 


17.8 ± 0.5% 


49.2 ± 1.0% 


77.0 ±2.1% 


85.8 ±5.1% 


1.0-2.0 


6.7 ± 0.5% 


7.9 ± 0.3% 


17.4 ±0.5% 


47.8 ± 1.1% 


76.2 ± 2.2% 


84.1 ± 5.6% 


2.0-100 


6.8 ±0.7% 


7.1 ±0.6% 


14.2 ± 0.7% 


37.2 ± 1.6% 


65.4 ± 3.4% 


83.1 ± 9.9% 



Note. — Breakdown of AGN selection efficiency by \og[LAGN/^-0 10 LQ] and E(B — V). Efficiencies were determined 
using Monte Carlo to generate model SEDs with varying contributions from the three normal galaxy templates, add AGN 
components to the model SED, and introduce photometric errors. Sec § 13. 2 1 for further details. 



